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Abstract. This paper extends parts of the results from |17l for plane curves 
to the case of hypersurfaces in R™. Let M be a compact connected oriented 
n — 1 dimensional manifold without boundary like S 2 or the torus S 1 X S 1 . 
Then shape space is either the manifold of submanifolds of R n of type M, or 
the orbifold of immersions from M to R n modulo the group of diffeomorphisms 
of M. We investigate almost local Riemannian metrics on shape space: These 
are induced by metrics of the following form on the space of immersions: 



where g is the standard metric on R n , f*g is the induced metric on M, h,k E 
C°° (M, R") are tangent vectors at / to the space of embeddings or immersions, 
where $ : R 2 — > R>o is a suitable smooth function, Vol(M) = J M vol( fg) 
is the total hypersurface volume of /(M), and where the trace Tr(L) of the 
Weingarten mapping is the mean curvature. For these metrics we compute 
the geodesic equations both on the space of immersions and on shape space, 
the conserved momenta arising from the obvious symmetries, and the sectional 
curvature. For special choices of <3> we give complete formulas for the sectional 
curvature. Numerical experiments illustrate the behavior of these metrics. 
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1. Introduction 

Many procedures in science, engineering, and medicine produce data in the form 
of shapes of point clouds in R n . If one expects such a cloud to follow roughly a 
submanifold of a certain type in R n , then it is of utmost importance to describe the 
space of all possible submanifolds of this type (we call it a shape space hereafter) 
and equip it with a significant metric which is able to distinguish special features 
of the shapes. Most of the metrics used today in data analysis and computer vision 
are of an ad-hoc and naive nature; one embeds shape space in some Hilbert space or 
Banach space and uses the distance therein. Shortest paths are then line segments, 
but they leave shape space quickly. 

Riemannian metrics on shape space itself are a better solution. They lead to 
geodesies, to curvature and diffusion. Eventually one also needs statistics on shape 
space like means of clustered subsets of data (called Karcher means on Riemannian 
manifolds) and standard deviations. Here curvature will play an essential role; 
statistics on Riemannian manifolds seems hopelessly underdeveloped just now. 

1.1. The shape spaces used in this paper. Thus, initially, by a shape we mean 
a smoothly embedded hypersurface in K. n which is diffeomorphic to M. The space 
of these shapes will be denoted B e = B e (M,M. n ) and viewed as the quotient (see 
[15] for more details in a more general situation) 

B e (M, R n ) = Emb(M,R n )/Diff(M) 

of the open subset Emb(Af, R n ) C C°°(M, R n ) of smooth embeddings of M in K", 
modulo the group of smooth diffeomorphisms of M. It is natural to consider all 
possible immersions as well as embeddings, and thus introduce the larger space 
Bi = Bi(M,R") as the quotient of the space of C°° immersions by the group of 
diffeomorphisms of M (which is, however, no longer a manifold, but an orbifold 
with finite isotropy groups, see [5]): 

Emb(M,K") »- Emb(M, R n ) / Diff(M) B e (M,R n ) 

p op 

III 
Imm(M, R n ) »- Imm(Af, R n )/ Diff(M) Bi(M,R n ) 

More generally, a shape will be an element of the Cauchy completion (i.e., the 
metric completion for the geodesic distance) of Bi(M, R n ) with respect to a suitably 
chosen Riemannian metric. This will allow for corners. On the other, discretizing 
for numerical algorithms will hide the need to go to the Cauchy completion. 

1.2. Where this paper comes from. In [17], Michor and Mumford have inves- 
tigated a variety of Riemannian metrics on the shape space 

Bi(S\R 2 ) = Imm(S ,1 ,M 2 )/Diff(S' 1 ) 

of immersion of the circle into the plane modulo the group of reparameterizations 
of the circle S 1 . In |T6J 3.10] they found that the simplest such metric has vanishing 
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geodesic distance; this is the metric induced by L 2 (arc length) on ImmfS 1 ,! 2 ): 

G° c (h,k)= f (h(9),k(6))\c'(6)\d0, 
Js 1 

c £ Imm(S' 1 ,R 2 ), h,ke C°°(S\R 2 ) = T c lmm(S\R 2 ). 

In [15] they found that the vanishing geodesic distance phenomenon for the L 2 - 
metric occurs also in the more general shape space Imm(M , N) / Diff (A/) where S 1 
is replaced by a compact manifold M and Euclidean K 2 is replaced by a Riemannian 
manifold N; it also occurs on the full diffeomorphism group Diff(iV), but not on 
the subgroup Diff (N, vol) of volume preserving diffeomorphisms, where the geodesic 
equation for the L 2 -metric is the Euler equation of an incompressible fluid. In |17[ 
section 3] a class of metrics was investigated which were called almost local metrics: 
they were of the form 

G?(M)= / Hl(c),K c {6))(h(e),k(6))\c'(6)\d6, 
Js 1 

where $ : M 2 — > R is a suitable smooth function, £(c) = f gl \c'(6) \ dO is the lenght 
of c, and k c is the curvature of c. If <f> = $(^(c)) then this is just a conformal 
change of the metric; it was proposed and investigated independently in |19) and 
in [331 (Ml [23 . For $ = 1 + Ak 2 the metric was investigated in great detail in [IB] . 

In this paper we take up the investigation of almost local metrics from [17] 
and we generalize it to the shape space Bi(M,R n ) = Imm(M, R")/ Diff(M) of 
hypersurfaces of type M in R ra ; here M is a compact orientable connected manifold 
of dimension n — 1, for example the hypersphere S n . 

1.3. Riemannian metrics. A Riemannian metric on Imm(M, K") is just a family 
of positive definite inner products Gf(h,k) where / <E Imm(A/, R n ) and h,k £ 
C°°(M,M. n ) = T f lmm(M 1 R n ) represent vector fields on R" along /. We require 
that our metrics will be invariant under the action of Diff (M), hence the quotient 
map dividing by this action will be a Riemannian submersion. This means that the 
tangent map of the quotient map Emb(Af, W 1 ) — > Bi(M,M. n ) is a metric quotient 
mapping between all tangent spaces. Thus we will get Riemannian metrics on Bi. 
Here the restriction to almost local Diff (M)-invariant metrics is very beneficial: For 
any / £ Emb(M, R n ) those vectors in T f Emb(M, IP) = C°° (M, R n ) which are G f - 
perpendicuar to the Diff (A/)-orbit through / are exactly those vector fields which 
are pointwise normal to f(M), i.e., T f {f o Diff (A/))- 1 ' ' = {h £ C°°(M,R n ) : 
h(x)±-T x f(T x M)}. We shall call such vectors horizontal. The tangent map of 
the quotient map Emb(Af, R") — > Bi(M,M. n ) is then an isometry when restricted 
to the horizontal spaces, just as in the finite dimensional situation. Riemannian 
submersions have a very nice effect on geodesies: the geodesies on the quotient 
space Bi are exactly the images of the horizontal geodesies on the top space Imm. 
(By a horizontal geodesic we mean a geodesic whose tangent lies in the horizontal 
bundle.) 

The simplest inner product on the tangent bundle to Imm(Af, R") is 
G° f (h,k)= [ g(h,k)vol(f*g), 
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where g is the Euclidean inner product on R n . Since the volume form vol(/*<?) 
reacts equivariantly to the action of the group Diff(M), this metric is invariant, 
and the map to the quotient Bi is a Riemannian submersion for this metric. 

All of the metrics we will look at will be of the form: 

G f (h,k)= f $(Yol(M),Tr(L))g(h,k)vol(fg) 

where $ : R 2 — s- R is a suitable positive smooth function, Vol(M) = J M vo\(f*g) is 
the total hypersurface volume of f(M), and where the trace Tr(L) of the Weingarten 
mapping is the the mean curvature. As mentioned above, these metrics will be 
called almost local. 



1.4. The contents of this paper. Let us now describe in some detail the contents 
of this paper and the metrics. 

Section [2] recalls concepts from differential geometry of hypersurfaces in a form 
that is suitable for our needs. 

In section [3] we calculate the derivatives of the metric, the volume form, the 
second fundamental form and some other curvature terms with respect to the im- 
mersion /, and the second derivatives for both tangent vectors horizontal. The 
differential calculus used is convenient calculus as in |llj . Some of the formulas can 
be found in [2 11 US], and in 0]. 

Section [4] describes the general Hamiltonian formalism that we shall use to com- 
pute geodesic equations and conserved quantities in a quick way. This is a shortened 
version of [T71 section 2], updated to the more general situation here. See also [T31 
section 2] for a detailed exposition in similar notation. First we consider general 
Riemannian metrics on the space of immersions which admit Christoffcl symbols. 
We express this as the existence of two kinds of gradients. Since the energy function 
is not even defined on the whole cotangent bundle of the tangent bundle we pull 
back to the tangent bundle the canonical symplectic structure on the cotangent 
bundle. Then we determine the Hamiltonian vector field mapping and, as a special 
case, the geodesic equation. We determine the equivariant moment mapping for 
several group actions on the space of immersions: The action of the reparametriza- 
tion group Diff(M), of the motion group of R", and also of the scaling group (if 
the metric is scale invariant). 



In section 5.1 we compute the geodesic equation on Imm(M, R") for the general 



almost local metric, and its conserved momenta. 



In section 6.2 we split this geodesic equation into its horizontal part (which is 
the geodesic equation on shape space Bi(M, R™)) and in its vertical part which 
vanishes identical. 

Section [7] computes the sectional curvature of Bi(M, R") for the general almost 
local metric. The result is the sum of terms P\P ... P$P and J M ~ . . . vol(^) 
over the terms QiQj for i = 1, . . . , 5, j = 1, . . . , 5. Some of these terms make pos- 



itive contribution to sectional curvature, some negative, and others are indefinite. 
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In section [8] some estimates for geodesic distance on Bi(M, R n ) are given: Two 
area swept out bounds, and Lipschitz continuity of ^/Vol(Af ). In addition we 
compare the almost local metrics to the Frechet metric. 

In section [9] we study the totally geodesic subspace of concentric spheres. 

In section [10] we specialize the general results to special choices for the weight 
function $. These are the following: 

• The G°-metric or L 2 -metric, where $ = 1, geodesic distance on Bi vanishes 
(as shown in |15j). and which has simple non- negative sectional curvature. 

• The G A -metric where $ = 1 + ATr(L) 2 ; its geodesic equation was treated 
in [16J for the situation of plane curves, with its sectional curvature. Here 
we have the geodesic equation and the sectional curvature as special cases. 

• Conformal metrics where $ = "I? (Vol). For curves these were investigated 
in [19j E«2 EH [25] . The full formula for sectional curvature is given only in 
the case that $(Vol) = Vol. 

• The scale invariant version where <& = Vol 1 -" +A v / . 

Section [Tl] contains numerical experiments for geodesies. We only do boundary 
value problems and no initial value problems (it is not clear that the equations 
are well posed). We use Mathematica to set up the triangulation of the surfaces, 
feed this into AMPL (a modelling software developed for optimization) and use the 
solver IPOPT. The numerical results are tested on the totally geodesic subspace of 
concentric spheres where we also have the analytic solutions. Then we study the 
translations of spheres for various metrics and discuss the appearing phenomena. 
Finally we deform some surfaces. 

This paper benefited from discussions with David Mumford, with Hermann 
Schichl who taught us about the use of AMPL, with Johannes Wallner and Tilak 
Ratnanather. 

Big parts of this paper can be found in the doctoral dissertation of Martin Bauer, 

2. Differential geometry of surfaces and notation 

In this section we will present and develop the differential geometric tools that 
are needed to deal with immersed surfaces. The most important point is a rigorous 
treatment of the covariant derivative and related concepts. 

In El section 2] one can find some parts of this section in a more general setting. 
We use the notation of [13] . Some of the definitions can also be found in [§J. A 
similar exposition in the same notation is [T]. 

2.1. Basic assumptions and convention. We always assume that M is a con- 
nected and compact manifold of dimension n — 1 without boundary. We will work 
with immersions of M into R™, i.e. smooth functions M — > M. n with injective 
tangent mapping at every point. We denote the set of all such immersions by 
Imm(M, K"). Immersions or paths of immersions are usually denoted by /. Vector 
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fields on Imm(M, E") or vector fields along / will be called h,k,m, for example. 
Subscripts like f t denote differentiation with respect to the indicated variable, so 
ft = dtf = df/dt, but subscripts are also used to indicate the foot point of a tensor 
field. 

2.2. Tensor bundles and tensor fields. We will deal with the tensor bundles 

T r s M T r s M ® f*TR n 

¥ 

M M 
Here TJM denotes the bundle of ( r s )-tensors on M, i.e. 

r s 

T r s M = (g)TM $ (g) T*M, 

and f*TR" is the pullback of the bundle TR n via /, see [H section 17.5]. A tensor 
field is a section of a tensor bundle. Generally, when E is a bundle, the space of its 
sections will be denoted by T(E). 

To clarify the notation that will be used later, some examples of tensor bundles 
and tensor fields are given now. 

• S k T*M = L* ym (TM;R) is the bundle of symmetric (I) -tensors, 

• A k T*M = L^ t (TM;R) is the bundle of alternating (l)-tensors, 

• n r (M) = T(A r T*M) is the space of differential forms, 

• X(M) = T(TM) is the space of vector fields, and 

• T(f*TR n ) = {he C°°(M, TR n ) : tt n o h = /} is the space of vector fields 
along f. 

2.3. Metric on tensor spaces. Let g E r(5?, T*K™) denote the Euclidean metric 
on R n . The metric induced on M by f € Imm(M, R n ) is the pullback metric 

g = f*g e T(S 2 >0 T*M), g(X, Y) = (f*g)(X, Y) = g(Tf.X, Tf.Y), 

where X, Y are vector fields on M. The dependence of g on the immersion / should 
be kept in mind. Let 

b = g : TM -> T*M and |) = cf x : T*M -> TM. 

g can be extended to the cotangent bundle T*M — T®M by setting 

g- 1 (a,[3) = g° 1 (a,(3) = a((3t) 

for a,P €T*M. The product metric 

r s 

extends g to all tensor spaces T r s M, and g r s ®g yields a metric on TJM ® f*TR n . 
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2.4. Traces. The trace contracts pairs of vectors and co-vectors in a tensor prod- 
uct: 

Tr : T*M ® TM = L(TM, TM) -»Mxl 

A special case of this is the operator ix inserting a vector X into a co-vector or 
into a covariant factor of a tensor product. The inverse of the metric g can be used 
to define a trace 

Tr 9 : T*M ® T*M ->Mxl 

contracting pairs of co-vecors. Note that Tr 9 depends on the metric whereas Tr 
does not. The following lemma will be useful in many calculations (see [H section 
2]): 

Lemma. 

g°(B, C) = Tt{g- x Bg- x C) for B,C £ T 2 °M i/B or C* is symmetric. 
(In the expression under the trace, B and C are seen maps TM — > T* M .) 

2.5. Volume density. Let Vol(M) be the density bundle over M, see [HI sec- 
tion 10.2]. The volume density on M induced by / <E Imm(M,E") is 

voi( ff ) = voi(/*§)er(Voi(M)). 

The volume of the immersion is given by 

Vol(/) = f vo\(f*g) = [ vol(g). 
Jm Jm 

The integral is well-defined since M is compact. If M is oriented we may identify 
the volume density with a differential form. 

2.6. Covariant derivative. We will use covariant derivatives on vector bundles 
as explained in [Ml sections 19.12, 22.9]. Let V 9 ,V 9 be the Levi-Civita covariant 
derivatives on (M, g) and (K™, g), respectively. For any manifold Q and vector field 
J on Q, one has 

:C°°(Q,TM) -> C°°{Q,TM), h^V 9 x h 

V| : C°°(Q,TR n ) -> C*°°(Q,T]R n ), h^V x h. 

Usually we will simply write V for all covariant derivatives. It should be kept in 
mind that V 9 depends on the metric g — f*g and therefore also on the immersion 
/. The following properties hold [HI section 22.9]: 

(1) Vx respects base points, i.e. n o Vj/i = tt o h, where 7r is the projection of 
the tangent space onto the base manifold. 

(2) Vx/i is C°°-linear in X. So for a tangent vector X x £ T X Q, Vx^h makes 
sense and equals iy x h){x). 

(3) Vxh is K-linear in h. 

(4) V x{o.h) = da(X).h + a.Vxh for a £ C° a (Q), the derivation property of 

(5) For any manifold Q and smooth mapping q : Q — > Q and Y y £ T y Q one 
has Vxq.Yyh = Vy (/ioq). If Y £ X(Qi) and X £ X(Q) are g-related, then 
Vy (h o q) = (Vxh) ° q- 
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The two covariant derivatives V x and V x can be combined to yield a covariant 
derivative Vx acting on C°°(Q, TJM®TR n ) by additionally requiring the following 
properties [H section 22.12]: 

(6) V x respects the spaces C°°(Q,T^M ® TN). 

(7) Vx(/i <£) k) = (Vx/i) ® k + h ® (Vjcfc), a derivation with respect to the 
tensor product. 

(8) Vx commutes with any kind of contraction (see [HI section 8.18]). A 
special case of this is 

V x («(y)) = (V x a){Y) +a(V x Y) for a ® Y : R n -> T^M. 

Property ([I]) is important because it implies that Vx respects spaces of sections of 
bundles. For example, for Q = M and / e C°°(M,R n ), one gets 

Vx : r(TJM ® fTR") -> r(TJM <g) f*TR n ). 

2.7. Swapping covariant derivatives. We will make repeated use of some for- 
mulas allowing to swap covariant derivatives. Let / be an immersion, h a vector 
field along / and X,Y vector fields on M. Since V is torsion-free, one has [HI 
section 22.10] 

(1) VxTf.Y- VyTf.X - Tf. [X, Y] = Tor(T/. A, Tf.Y) = 0. 
Furthermore one has [HI section 24.5] 

(2) VxW/i - VyVxft - V [x , Y ]h = R" ° (Tf.X, Tf.Y)h = 0. 

These formulas also hold when / : R X M — > R" is a path of immersions, 
h : R x A/ — > TM" is a vector field along / and the vector fields are vector fields on 
R x M. A case of special importance is when one of the vector fields is (dt, 0m) and 
the other (Ok, Y), where Y is a vector field on M. Since the Lie bracket of these 
vector fields vanishes, ([!]) and ^ yield 

(3) V (at ,o M )T/.(0 R ,y) - V (OK , Y) Tf.(d t ,0 M ) = 
and 

( 4 ) v (9 t ,OM)V ( o E ,i-)/i - V(o K ,r)Vj;a tj o M )^ = 0. 

2.8. Higher covariant derivatives and the Laplace operator. When the co- 
variant derivative is seen as a mapping 

V : r(TJM) -> r(IJ +x M) or V : T{T r s M ® f*TR n ) -> T(T r s+1 M ® f*TR n ), 

then the second covariant derivative is simply VV = V 2 . Since the covariant 
derivative commutes with contractions, V 2 can be expressed as 

V| Y := l y l x V 2 = tyVxV = VxW - V Vx y for X, Y € X(M). 

Higher covariant derivates are defined as V fc , k > 0. We can use the second co- 
variant derivative to define the Laplace- Bochner operator. It can act on all tensor 
fields B and is defined as 

AB = -Tr 9 (V 2 B). 
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2.9. Normal bundle. The normal bundle Nor(/) of an immersion / is a sub- 
bundle of f*TR n whose fibers consist of all vectors that are orthogonal to the 
image of /: 

Nor(/) x = {Y G T f(x) N : VX G T X M : g(Y,Tf.X) = 0}. 

Any vector field h along / can be decomposed uniquely into parts tangential and 
normal to / as 

h = Tf.h T + h ± , 

where h T is a vector field on M and h 1 - is a section of the normal bundle Nor(/). 
When / is orientable, then the unit normal field v of / can be defined. It is a 
section of the normal bundle in one of the above forms with constant (/-length one 
which is chosen such that 

{v{x), T x f.X u T x f.X 2 , . . . T x f.X n ^) 

is a positive oriented basis in Tf/ x \M. n if X\,. . . ,X n _i is a positive oriented basis 
in T X M. In this notation the decomposition of a vector field h along / reads as 

h = Tf.h T + a.v. 

The two parts are defined by the relations 

a = g{h,v) G C°°(M) 

h T G X(M), such that g{h T ,X) = g(h,Tf(t,-).X) for all X e X(M). 

2.10. Second fundamental form and Weingarten mapping. Let X and Y be 

vector fields on M. Then the covariant derivative VxTf.Y splits into tangential 
and a normal parts as 

VxTf.Y = Tf.(V x Tf.Y) T + (VxTf.Y^ = T.f.V x Y + S(X,Y). 

S is the second fundamental form of f. It is a symmetric bilinear form with values 
in the normal bundle of /. When Tf is seen as a section of T* M ® f*TR n one has 
S = VTf since 

S(X, Y) = V x Tf.Y - Tf.V x Y = (WTf)(X, Y). 

Taking the trace of S yields the vector valued mean curvature 

Tr s (5) G r(Nor(/)). 

One can define the scalar second fundamental form s as 

s(X,Y) = g(S(X,Y),u). 

Moreover, there is the Weingarten mapping or shape operator L = g~ l s. It is a 
(/-symmetric bundle mapping defined by 

s(X,Y)=g(LX,Y). 

The eigenvalues of L are called principal curvatures and the eigenvectors principal 
curvature directions. Tr(L) = Tr s (s) is the scalar mean curvature and for surfaces 
in K 3 the Gaufi- curvature is given by det(L). The covariant derivative V xv of the 
normal vector is related to L by the Weingarten equation 

V x v = -Tf.L.X. 
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2.11. Directional derivatives of functions. We will use the following ways to 
denote directional derivatives of functions, in particular in infinite dimensions. 
Given a function F(x, y) for instance, we will write: 

D( x .h)F as shorthand for dt\oF(x + th,y). 

Here (x, h) in the subscript denotes the tangent vector with foot point x and di- 
rection h. If F takes values in some linear space, we will identify this linear space 
and its tangent space. 



3. Formulas for first variations 

Recall that many operators like 

g = f*g, S = Sf, vol(g), V = V», A = A 9 , ... 

implicitly depend on the immersion /. We want to calculate their derivative with 
respect to /, which we call the first variation. We will use this formulas to calculate 
the metric gradients that are needed for the geodesic equation. 

Some of the formulas can also be found in (2J SI HSl EU H] ■ 

3.1. Paths of immersions. All of the concepts introduced in section [2] can be 
recast for a path of immersions instead of a fixed immersion. This allows to study 
variations immersions. So let / : K — » Imm(M, N) be a path of immersions. By 
convenient calculus llj, / can equivalently be seen as / : K. X M — > N such that 
f(t, •) is an immersion for each t. We can replace bundles over M by bundles over 
IxM: 

pr^ TJM pr^ TJM ® f*TN Nor(/) 

T ¥ Y 

Rx M Rx M Rx M 

Here pr 2 denotes the projection pr 2 :MxM-> M. The covariant derivative V zh is 
now defined for vector fields Z on R x M and sections h of the above bundles. The 
vector fields (9(,0m) and (Or,X), where X is a vector field on M, are of special 
importance. Let 

ins t :M->Rx M, x ^ (t, x). 
Then by Q3 22.9.6] one has for vector fields X, Y on M 

V x Tf{t, -).Y = V x T{f o ins t ) o Y = V x Tf o Tins t oY 

= V X Tf o (Or, Y) o ins t - V T ins t oxTf o (Or, Y) 

= (V ( o K ,x)T/o(0 K ,y))oins t . 

This shows that one can recover the static situation at t by using vector fields on 
R x M with vanishing M-component and evaluating at t. 
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3.2. Setting for first variations. In all of this chapter, let / be an immersion 
and ft € Tt Imm a tangent vector to /. The reason for calling the tangent vector 
ft is that in calculations it will often be the derivative of a curve of immersions 
through /. Using the same symbol / for the fixed immersion and for the path of 
immersions through it, one has in fact that 



For the sake of brevity we will write dt instead of (94,0m) and X instead of (Or, X), 
where X is a vector field on M. 

Let the smooth mapping F : Imm(M, N) — > T(TJM) take values in some space 
of tensor fields over M, or more generally in any natural bundle over M, see [ID] . 

3.3. Lemma (Tangential variation of equivariant tensor fields). If F is equivariant 
with respect to pullbacks by diffeomorphisms of M, i.e. 



for all ip £ Diff(Af) and f € Imm(M, N), then the tangential variation of F is its 
Lie-derivative: 



This allows us to calculate the tangential variation of the pullback metric and 
the volume density, because these tensor fields are natural with respect to pullbacks 
by diffeomorphisms. 

3.4. Lemma (Variation of the metric). The differential of the pullback metric 



Proof. Let / : R x M — > N be a path of immersions. Swapping covariant derivatives 
as in section |2.7| formula pi) one gets 



d t (g(X,Y)) = d t (g(Tf.X,Tf.Y)) - g(V dt Tf.X,Tf.Y)+g(Tf.X,V dt Tf.Y) 
= g(Vxft,Tf.Y)+g(Tf.X,V Y ft) = (2 Symg(V/t, Tf))(X, Y). 



D UJt) F = d t F(f(tj). 



F(f) = (ip*F)(f) = <p*(F((<p- 1 yf)) 



D (f,Tf.fT ) F = d t \ Q F(foFlf) =d t \ Q F((Flf )*/) 
= 8 t \o(Flf)\F(f))=C f T(F(f)). 



Imm -> r(S* T*M), 
f >-> 9 = f*9 



is given by 



D (fJ t )9 = -2g(Jt,v).s + £ f r {g). 




D (f-Tf.fJ)9 = 2SymV(/ t T ) b = C f rg 



follows from the equivariance of g (see 3.3). 



□ 
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3.5. Lemma (Variation of the inverse of the metric). The differential of the inverse 
of the pullback metric 

Imm -> T{L{T*M,TM)), 
is given by 

Duj^g- 1 = -2ff(/ t ,^).L. ff - 1 +£ / T( 9 - 1 ). 

Proof. 

dtg- 1 = -g'Hd^g- 1 = -g~ l {- 2g{ft,S)+C f7 g)g- 1 

= 2g- 1 g(f^ S)^ 1 - g- 1 (C f rg)g- 1 = 2g(f t ± ,g- 1 Sg- 1 ) + ^(g- 1 ) □ 

3.6. Lemma (Variation of the volume density). The differential of the volume 
density 

Imm ->■ Vol(M), 
/ ^ vo%)=vol(/*s) 

is given by 

D {fJt) vol(g) = ( div*(/ t T ) - s(/ t V). Tr(L)) vol(g). 

Proof. Let g(t) e r(5'?,gT*M ) be any curve of Ricmannian metrics. Then 

d t vol(g) = ^Tr(g- 1 .d t g)vol(g). 
This follows from the formula for vol(g) in a local oriented chart (u 1 , . . . u n ) on M: 
d t vol(g) = a t y 1 Aet{{gij) l3 ) du 1 A • • • A du™" 1 

= ~ Ta = , , Tr(adj(. 9 )a tff ) du 1 A • • • A d^ 1 

' Tr{det{{g l0 ) lJ )g- 1 d t g) du 1 A • • • A du"- 1 



2- v /det((<fa)^, 
= l -Tr{g- 1 .d t g)vo\(g) 

Now we can set g — f*g and plug in the formula for d t g = d t (f*g). This yields 
S t vol(<7) = \Tv{g-\-2-g{f uV ).s + C h T 9 )).yo\{g) 

= -gifu^Trig-Ks^volig) + \ Tr^" 1 ^). vol(g) 
The same calculation as above with d t replaced by £ h r shows that 
C h r vol(g) = \ Tx^L^g). vo%). 

Therefore 

t vol(g) = -g(f t ,v) Tr(L). vol(g) + C h r (vol(g)) 

= -g{f u v) Tr(L). vol(<?) + div^ T ) vol( 5 ). □ 
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3.7. Lemma (Variation of the volume). The differential of the total Volume 

Imm — > R, 
/ ^ Vol(f) = f M vol(f*g) 



is given by 



D(f Jt ) vol(g) = - / g{ft, u). Tr(L) vol(«?). 

JM 



3.8. Lemma (Variation of the second fundamental form). The differential of the 
second fundamental form 

f Imm T(S 2 T*M), 
{ f -> at 

is given by 

D (fM s = .9(V 2 / t ,^) = V 2 .g(/ t , v) - g(f t , v).g o (L <g> L) + £ /t T. S 

Proo/. By definition s(X,Y) = g(S(X,Y),u) = gjV x (Tf.Y) - Tf.V x Y, v) . Inter- 
changing covariant derivatives as in pTflpl) andpTrlflil yields: 

d t s(X, Y) = g{d t S{X, Y),u) + g(S(X, Y),d t v) 

= giV x V Y Tf.dt - V VxY Tf.d t , v)+0 = g(V 2 xx f t , v) 

where the term g{S{X, Y), dtv) vanishes since dtv is tangential (see 3.11 1. For the 
normal part this yields: To get the second formula we calculate: 

D (f,9(fw)-v) s ( X > Y ) = 9{^x,y{9Uuv).v),v} 

= ^x,y9{M ~ Sift, v).g^ x v, V Y v) + 

= V^ y s(/ t) i/) - gift, y).g{LX, LY) 

By |3.3[ the formula for the tangential variation follows from the equivariance of the 
second fundamental form: 

sf*(X, Y) = g(V x T(/ o <j>) o Y, v^) = g(V x (Tf o (<f> m Y) o <f>), v* o cj>) 

= g(^T4>oxTf o (<p*Y),l/f o (f) o (f) 

= g(W$,xTf o i^Y), vf)o(f > = affaX, </>*Y) o = {4>*a f ){X, V) □ 

3.9. Lemma (Variation of the Weingarten map). XTie differential of the Weingarten 
map 

Imm -> r(End(TM)), 



f ^ Lf 

is given by 

D (f Jt) L = g-\V 2 {gift, g(f u v)L 2 + £ f riL). 

Proof. From L = g .s follows 
d t L = .g- 1 .d iS + 9 t ( ff - 1 ). S 

= g-'.^igifuv)) - gift,v)g.L 2 + £ f ris)) + faft^Lg- 1 +£ / t(<t 1 )). 
= g-W&ft, v)) + g(fu y)-L 2 + % (L). □ 
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3.10. Lemma (Variation of the mean curvature). The differential of the mean 
curvature 

Imm -> C°°(M), 
f » Tr(Lf) 
is given by 

D (fJt) Tr(L) = -A(g(f t ,v)) + g(f u v).Tr (L 2 ) + d( Tr(L)) (/ t T ). 

Proof. This statement follows from the linearity of the trace operator and from the 
previous equation for D(jj t )L. □ 

3.11. Lemma (Variation of the normal vector field). The normal vector field is a 
smooth map v : R x M — > R n . Therefore, as explained in section \2.6[ we can take 
its covariant derivative along vector fields onix M. Identifying d t with the vector 
field (dt, 0m) on R x M, we get 

V 9t ^ = -T/.(L/ t T +grad 9 (g(f t ,v)j). 

Proof. V ' Q t v is tangential because g(VQ t v,v) = \dtg(y,v) = 0. Therefore one can 
write V dt v = Tf.(y dt v) T . Then for all X € X(M) we have 

g{{V dt v) T ,X)=g(V dtV ,Tf.X) = - g(v,V dt TfX) = -g{v,V x Tf.d t ), 

where in the last step we swapped X and dt as in section 2.7 formula (J3|. Splitting 
into normal and tangential parts yield: 

g((V 9t u) T ,X) = -g{u,Vxh) = -§(?, Vx{Tf.fJ + g{ft,v).v)) 

= -g(v,V x (Tf.f? +g{fuv).vj) 
= -s{X,fJ)-V x {g{f t ,v))-0 

= -g(LfJ +gr A tfg{ft,v),x) □ 

In this section, let V = V 9 = V-^ 5 be the Levi-Civita covariant derivative 
acting on vector fields on M. Since any two covariant derivatives on M differ by a 
tensor field, the first variation of V-^ s is tensorial. It is given by the tensor field 
D(f,f t )V f ' s G r(T'M). 

3.12. Lemma (Variation of the covariant derivative). The tensor field -D(/,/ t )V-^ s 
is determined by the following relation holding for vector fields X,Y, Z on M: 

g({D UJt) \7)(X,Y),Z) = ~(yD(f ift) g)(X®Y®Z + Y®X®Z-Z®X®Y) 

Proof. The defining formula for the covariant derivative is 

g{V x Y, Z)= l - [X 9 (Y, Z) + Yg(Z, X) - Zg(X, Y) 

- g(X, [Y, Z\) + g(Y, [Z, X]) + g(Z, [X, Y}) 
Taking the derivative D(fj t ) yields 

(D UJt) g)^xY, Z) + g((D {fJt) V)(X, Y), Z) 
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1 r X((D Utft) g)(Y,Z)) +Y({D UJt) g)(Z,X)) - Z((D {fJt) g)(X,Y)) 
- (D {fJt) g)(X, [Y,Z]) + (D UJt) g)(Y, [Z,X]) + (D (fJt) g)(Z, [X,Y]) 

Then the result follows by replacing all Lie brackets in the above formula by co- 
variant derivatives using [X, Y] — Vx7 — VyX and by expanding all terms of the 
ioxm X{{D UJt g){Y,Z)) using 

X((D(fj t )9)(Y,Z)) = 

(\7 x D (fJt) g)(Y,Z) + (D (ftft) g)(VxY,Z) + (D {fJt) g)(Y,\7 x Z). □ 

3.13. Setting for second variations. All formulas for second derivatives will be 
used in section \7. 21 There we consider a curve of immersions 

f(t,x) = exp Mx) {t.a(x).v fa (x)) = f (x) + t.a(x).v h (x), 

for a fixed immersion /o. This curve of immersions has the property that at t = 
its first derivative and the covariant derivative of the first derivative are both 
horizontal, i.e., 

(1) /| t =o = /o, d t \vf = a.v f \ and V dt Tf.d t \ t=Q = 0. 

In all calculations of second variations we will assume that the above properties 
hold. 

3.14. Lemma (Second variation of the metric). The second derivative of the pull- 
back metric 

Imm T(Sl T*M), 
f -> 9 = f*9 



along a curve of immersions f satisfying property ([!]) from section 3.13 is given by 
df\ f*g = 2{da ® da) + 2a 2 g o (L /o ® L fo ). 

Proof. Since Vd t T/.<9 t | = 0, we have 

d?\ g(X,Y)=d?\ g(Tf.X,Tf.Y) 

= d t \ g(V 9t Tf.X,Tf.Y)+d t \ g(Tf.X,V dt Tf.Y) 
= 2g(V 9t Tf.X\ , V dt Tf.Y\ ) + + = 2g(V x Tf.d u V Y Tf.d t ) 
Using Tf.dt = a.v*° we get 

d?\ (g(X,Y)) = 2da(X).da(Y) + 2a 2 g{V x v h ,V Y v fa ) 

= 2(da®da)(X,Y) + 2a 2 .g (L fo X, L fo Y) □ 

3.15. Lemma (Second variation of the inverse metric). The second derivative of 
the inverse of the pullback metric 

' Imm -+ T(L(T*M,TM)), 



along a curve of immersions f satisfying property ([!]) from section 3.13 is given by 
dMfgr 1 = Qa 2 (Lfo) 2 .g^ - 2g^ (da ® da)g^ . 
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Proof. We look at g = f*g as a bundle map from TM to T*M. Then 

#t lo(s _1 ) = dt\ a (- g^Ag-g^ 1 ) = 2g^ .d t \og -g^ 1 -d t \og -g^ 1 - sb 1 - 9 ? los-ffo 1 

= 2{-2aL fa ) 2 .g^ 1 - g^ 1 .(2(da ® da) + 2a 2 g o (L /o ® L 1 "))^ 1 

= 8a 2 (L fo ) 2 .go 1 - 2g^ 1 (da ® dajg^ 1 - 2a 2 (L f ") 2 .g^ 1 □ 

3.16. Lemma (Second variation of the volume form). The second derivative of the 
volume form 

( Imm -> fT^M), 
\ / >-> vol(g) = vol(/*3) 
along a curve of immersions f satisfying property ([!]) from section 3.13 is given by 

d 2 | vol(. 9 ) = \a 2 Ti(L^ ) 2 — a 2 Tr {[L^ ) 2 ) + \\da\\ 2 } vo\(g ), 



Proof. In section 3.6 we showed that for any curve of Riemannian metrics g(t) <E 



r(S 2 T*M), we have 



Therefore 



^vo%) = -Tr(.g- 1 .a t5 )vol(.g). 



d 2 vol(g) = d t \ Tr {g' 1 .d t g) vol(g) = - Tr {d t (g- x ) .d t g) vol(.g) 

+ - Tr (g-Ktfg) vol(g) + ~ Tr (g-Kdtg)* vol(g) 
Evaluating at f = and setting g(t) = f*g we get 

d 2 \ vo\(g) = l - Tr {{2aL^ g^).(^2a.s^)) vol( 5o ) 
+ - Tr (gQ 1 .2(da ® da)) vol(g ) 
+ iTr(.g - 1 .2a 2 ffo .(L^) 2 )vol( 5 o) 
+ 1 Tr ( 5o - 1 .(-2a.s^))(- Tr(i^).a) vol(«? ) 
a 2 Tr(L'°) 2 - a 2 Tr ((i>) 2 ) + \\da\\% } vo\(g ) 



a 



3.17. Lemma (Second variation of the second fundamental form). The second 
derivative of the second fundamental form 

Imm T(S 2 T*M), 
f * J 

along a curve of immersions f satisfying Property \3. 13$ 1 ) is given by 



d 2 \ s = 2{da® da)(ld®L fo + L fo ® Id) - ||cia||^-i .s fo 



O t s(X, Y) = g(y 2 XY Tf.d u v) = g(V x<Y f u v). 



+ 2.a(V grad3 o( a )S /o ). 
Proof. From section [378] we have 
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Using Va t /* = we get 

dfs(X,Y) = g{V 2 xy f u V dt y)+g{V dt V x V Y f t ~V dt V VxY f u y) 
= 9{^ 2 x.Yh^d t v) + o - 9(Vv x YV a Jt + V [dt}VxY] f u v) 
= 9(Vx,r/t, Vg,i/) + - g(V [ai , Vx y]/ f , i/) 

ifJt )V)(X.Y)ft, V) 

In the last step we used 

[d t ,vpy] = p t ,o M ),(o R ,vpy)] 

= V)(x,y)) - (d (Wi) v)(it). 

Evaluating at £ = yields: 

a t 2 | oS (x,r) = 

= ff(V^jr(0.I/A), _ T/q . grad 90 fl) _ j(V (J3(ft) ^ ft)V)(Xfl (fl>),^) 

= + .g(rfa(X).Vyi/ /o +da(r).Vx^ /o ,-T/ .grad S0 a) 
+ ff(a.V^(i>), -T/ . grad 90 a) - da{(D {fo ^ fo) V)(X, V)) + 0. 
We will treat the three terms separately. The first one gives, using V z v = —Tf.L.Z: 
g(da{X).V Y is fo + da(Y)S7 x is f « , -T/ . grad 90 a) = 
= .9o(dapO£ /o ^ + da(Y)L fo X : grad 90 a) 
= da(X).da(L fo Y) + da(Y).da(L fo X). 
For the second term we get: 
ff(a.V^ y (i/A), -T/ . grad 90 a) = -ag{V x V Y v la - V Vx >V , T/ . grad 90 a) 
= -ag(-W x (Tf L f °Y) + Tf L f "V x Y, Tf . grad 90 a) 
= -ag(-(VTf )(X, L So Y) - Tf a V x (L f °Y) + TfoL fa V x Y, Tf . grad 90 a) 
= + ag(Tf (V x L fo )(Y),Tf . grad 90 a) = a 5o ((V x £ /o )(n g^d 90 a) 
= aVx(9o(£ /o y,grad s ° a)) - a^i/ V^F, grad 90 a) 

-a. 9o (i /o ^V x grad 90 a) 
= aV x (s /o (F, grad 90 a)) - as /o (V x 7, grad 90 a) - as /o (Y, V x grad 90 a) 
= a(V x s)(T,grad 90 a) 

2 

V x y ^ is symmetric in X,Y because the ambient space R™ is flat. Therefore the 
last formula and the symmetry of s imply that 

a(Vxs)(y,grad 90 a) = o(V y s)(X, grad 90 a) = a(V grad90 a s)(X, Y). 

The third term yields, using formula |3.12| 

-.9o(P (/ , a .,/o)V)(X,y),grad 9o (a)) = 

= -i(V(-2a.s /o ))(X,y,grad 9o (a)) - i(V(-2a.s /o ))(F, X, grad 9o (a)) 
+ ^(V(-2a. S /°))(grad 90 (a),X,y) 
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= da{X).s fo (Y,gr&d 9o (a)) + a.(V x s fo )(Y, gr&d 90 (a)) 

+ da(Y).s fo (X, grad 90 (a)) + a.(Vy s fo ) {X, grad 90 (a)) 

- da(grad 9o (a))>(X, F) - «-(V gradS o( Q )S /o )(X, F) 
= da(X).da(L fo Y) + da(Y).da(L fo X) 

- \\daW]., .s fo (X,Y) + a.(V grad9 o( a)S /o )(X,y). □ 

3.18. Lemma (Second variation of the mean curvature). The second derivative of 
the mean curvature 

( Imm -> C°°(M), 
\ / -> Tr(L/) 
along a curve of immersions f satisfying Property f3.l3[ l ) is given by 

<9 t 2 | Tr(L) = 2a 2 Tr((L /o ) 3 ) + AaTi (L fo g^ 1 .V 2 a) + 2Tr(g- 1 (da® da)L fo ) 
- \\daf g -i Tr(L fo ) + 2aTr 9 ° (V grad30 a s fo ) 

Proof. From Tr(L) = T¥(g .s) we get 

9 2 Tr(L) =Tr(a?(o- 1 ). S ) + 2Tr (c^g" 1 ).^) + Tr (g" 1 .^*) 



Evaluating at t — we get 

<9 t 2 |o Tr(L) = Tr (6a 2 [L fo f .g^ 1 J°) + Tr ( - 2. 9( 7 1 .(da ® do).flo X .« A ) 
+ 2Tr (2aL f "g^ 1 .\7 2 a) + 2Tr {2aL h g^ 1 .(-ag Q {L fo f)) 
+ 2. Tr (flfo ^((da ®dao i>) + (da o L /o fg) da)) 

- ||da||J-i Tr(Z/°) + 2a Tr 90 (V gradso a s fo ) 
= 2a 2 Tr ((£ /o ) 3 ) - 2Tr (^.(da ® da).L fa ) 

+ 4a Tr (L /o 1 . V 2 a) + 4 Tr (.a^ 1 . (da <g> da) . (L /o )) 

- ||da||J-i Tr(Z/°) + 2a Tr 90 (V gradso a s fo ) □ 

4. The Hamiltonian approach 

See |13[ section 2] for a detailed exposition in similar notation. 

4.1. The setting. Consider the smooth Frechet manifold Imm(M, R") of all im- 
mersions M — > M. n , which is an open subset of C°°(M, M. n ). 

The tangent bundle is TTmm(M, R n ) = Imm(M,R n ) x C°°(M,R n ), and the 
cotangent bundle is T* Imm(M,R") = Imm(M,R n ) x V(M) n where the second 
factor consists of n-tuples of distributions in T>(M) — C°°(M)' which is the space 
of distributional sections of the density bundle. 

We consider smooth Riemannian metrics on Imm(Af, R n ), i.e., smooth mappings 
G : Imm(M,R™) x C°°(M,R n ) x C°°(M,R n ) R 
(/, h, k) i — y Gf(h,k), bilinear in h, k 
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G f (h,h)>0 for h 7^0. 

Each such metric is weak in the sense that G/, viewed as bounded linear mapping 

Gf : T f Imm(M, R n ) = C°°(M, R") -> T} Imm(M, R") = V(M) n 

G : TImm(M, R n ) -> T* Imm(M, R") 

G(f,h) = (f,G f (h, .)) 

is injective, but can never be surjective. We shall need also its tangent mapping 

TG : T(rimm(M,R n )) ->■ T(T* Imm(M,M n )) 

We write a tangent vector to TImm(M, R n ) in the form (/, h; k, V) where (/, h) e 
T Imra(M, M. n ) is its foot point, k is its vector component in the Imm(M, Re- 
direction and where V is its component in the G°°(M, R™)-direction. 

Then TG is given by 

TG(f,h;k,V) = (f,G f (h, .);k,D (m G f (h, .) + G f (V, .)) 

Moreover, if X = (/, h; k, V) then we will write X\ = k for its first vector component 
and X 2 = V for the second vector component. 

Note that only these smooth functions on Imm(M,R n ) whose derivative lies in 
the image of G in the cotangent bundle have G-gradients. This requirement has 
only to be satisfied for the first derivative, for the higher ones it follows (see (llj). 
We shall denote by Cq (lmm(Al, R")) the space of such smooth functions. 

We shall always assume that G is invariant under the reparametrization group 
Diff(M), hence each such metric induces a Riemann- metric on the quotient space 
Bi(M, R") = Imm(M,R")/Diff(M). 

In the sequel we shall further assume that that the weak Riemannian metric G 
itself admits G-gradients with respect to the variable f in the following sense: 



D (ftm) G f (h, k) = G f (m, H f (h, k)) - G f (K f (m, h),k) 



where 



H,K : Imm(M,R n ) x G°°(Af,R") x G°°(M,R n ) -> G°°(M,R n ) 

(f,h,k) i ^ Hf{h,k),K f (h,k) 

are smooth and bilinear in h,k. 

Note that H and K could be expressed in abstract index notation as gij,kg kl and 
gij,kg zl - We will check and compute these gradients for several concrete metrics 
below. 

4.2. The fundamental symplectic form on TImm(M,R n ) induced by a 
weak Riemannian metric. The basis of Hamiltonian theory is the natural 1- 
form on the cotangent bundle T* Imm(M, R") given by: 

9 : T(T* Imm(M,R™)) = Imm(M,R") x V(M) n x C°°(M, R") x V(M) n -> R 

(f,a;h,P) (a,h). 

The pullback via the mapping G : TImm(M, R") -> T* Imm(Af, R n ) of the 1-form 
Q is then: 

(G*Q) {fth) (f,h;k,V) = G f (h,k). 



20 



MARTIN BAUER, PHILIPP HARMS, PETER W. MICHOR 



Thus the symplectic form uo = —dG*® on TImm(M, M. n ) can be computed as 
follows, where we use the constant vector fields (/, h) h- > (/, h; k, V): 

W(/,/ l )((fci,y 1 ),(fc 2 ,F 2 )) = -d(G*6)((fc 1 ,y 1 ),(A :2 ,V r 2 ))| (/i , i) 

= -D (fM) G f (h, k 2 ) - G S {V X M) + D UM) G f (h, h) + G f {V 2 ,k x ) 

(1) = G f (k 2 ,H f (h, h) - K f (k u h)) + Gf(V 2 , h) - G f {V u k 2 ) 

4.3. The Hamiltonian vector field mapping. Here we compute the Hamilton- 
ian vector field gra,d LU (F) associated to a smooth function F on the tangent space 
TImm(A/,r), that is F € Gg?(Imm(Af, R n ) x C°°(M,R n )) assuming that it has 
smooth G-gradients in both factors. See [HJ section 48]. 

Using the explicit formulas in |4.2[ we have: 
cj (m (grad w (F)(/ ; h),(k,V)) 

= u, (fth) ((grad? h),grad%(F)(f, hj), (k, V)) = 
= G f (k,H f (h, grad? (F) (f, h)))-G f (K f (grad? (F) (f,h),h),k) 
+ Gf (V, grad? (F) (/, h))-G f (grad^ (F) (f,h), k) 
On the other hand, by the definition of the w-gradient we have 
co (W (grad"(F)(/, h), (k, V)) = dF(f, h)(k, V) = D (m F(f, h) + D {hjV) F(f, h) 
= G f (gradf (F) (f,h),k) + Gf (grad 2 G (F) (/, ft) , V) 
and we get the expression of the Hamiltonian vector field: 



grad^(F) (f,h) = 


gradf (F)(/» 


grad^(F)(f,h) = 


-gradf (F)(f, h) 




+ H f (h, gradf (F)(/, h)) - (gradf (F)(/, h),h) 



Note that for a smooth function F on TImm(M, K n ) the w-gradient exists if and 
only if both G-gradients exist. 

4.4. The geodesic equation. The geodesic flow is defined by a vector field on 
TImm(Af, M"). One way to define this vector field is as the Hamiltonian vector 
field of the energy function 

E(f,h) = -G f (h,h), E : Imm(Af,E n ) x G°°(Ai",R n ) ->■ E. 
The two partial G-gradients are: 

G/ (gradf (£)(/, h),V) = d 2 E(f, h){V) = G f (h, V) 

gradf(£)(/,/i) = /i 
G f (gradf (£?)(/, h),k) = d x E(f, h)(k) = ±D (m G f (h, h) 
= \Gf{k,H f {h,h)) 
gr a df(E)(f,h) = ±H f (h,h). 
Thus the geodesic vector field is 

grad? (£)(/, h) = h 
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grad£ (£)(/, h) = \H f {h, h) - K f (h, h) 
and the geodesic equation becomes: 

\h t ZlH f (h,h)-K f (h,h) ° r \ f" = hHf{ft,ft)-K f {f t ,f t ) 

This is nothing but the usual formula for the geodesic flow using the Christoffel 
symbols expanded out using the first derivatives of the metric tensor. 

4.5. The momentum mapping for a G-isometric group action. We consider 
now a (possibly infinite dimensional regular) Lie group with Lie algebra g with a 
right action g i— > r 9 by isometries on Imm(M, R"). Denote by X(Imm(M, W" 1 )) the 
set of vector fields on Imm(M , E"). Then we can specify this action by the funda- 
mental vector field mapping Q : g —> £(Imm(M,R n )), which will be a bounded Lie 
algebra homomorphism. The fundamental vector field Qx , X € g is the infinitesimal 
action in the sense: 

Cx(f) = d t \ a r^ tx \f). 
We also consider the tangent prolongation of this action on TImm(M, R n ) where 
the fundamental vector field is given by 

CF mm : (f,h) ^ (f,h;Cx(fhD (w (Cx)(f) =: &(/,/>)) 
The basic assumption is that the action is by isometrics, 

G f (h,k) = ((r')*G) / (M) = G rS(/) (T / (r»)/ ll T J (r»W. 
Differentiating this equation at g — e in the direction X £ g we get 

(1) = D (Mx(f)) G f (h,k) + G f (C x (f,h),k) + G f (h,Cx(f,k)) 

The key to the Hamiltonian approach is to define the group action by Hamiltonian 
flows. We define the momentum map j ; g — > Cg 3 (TImm(M, W 1 ), R) by: 

jgZ h) = G f (t x (fihy 

Equivalently, since this map is linear, it is often written as a map 

J :Tlmm(M,R n )^g', (J(f,h),X)=j x (f,h). 

The main property of the momentum map is that it fits into the following commuta- 
tive diagram and is a homomorphism of Lie algebras: 

H (T Imm) ^ C% (T Imm, R) > X(T Imm, w) H 1 (T Imm) 






where I(TImm,w) is the space of vector fields on TImm whose flow leaves u> fixed. 
Note also that J is equivariant for the group action. See [T7] for more details. 

By Emmy Noether's theorem, along any geodesic t t— > c(t, . ) this momentum 
mapping is constant, thus for any X £ q we have 

(J(f,ft),X)=j x (J,ft) = G f (Cx(f),ft) is constant in t. 



We can apply this construction to the following group actions: 
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• The smooth right action of the group Diff(M) on Imm(M,R n ), given by 
composition from the right: / h-> / o ip for ip £ Diff (M). 

For X £ X(M) the fundamental vector field is then given by 

Cx m (/) - Cx(f) = x(f) = dMf ° Fif ) = df.x 

where Fl t denotes the flow of X. The reparametrization momentum, for 
any vector field X on M is thus: 

j x (f,h) = G f (df.X,h). 

Assuming the metric is reparametrization invariant, it follows that on any 
geodesic f(x,t), the expression Gf(df.X,f t ) is constant for all X. 

• The left action of the Euclidean motion group E" x SO(n) on Imm(M,M n ) 
given by / i-> A/ + B for (-B, A) € R" x SO{n). The fundamental vector 
field mapping is 

C(B,x)(f) = Xf + B 

The linear momentum is thus Gf(B, h), B 6 R n and if the metric is transla- 
tion invariant, Gf(B, f t ) will be constant along geodesies for every B £ M. n . 
The angular momentum is similarly Gf(X.f, h), X £ so(n) and if the metric 
is rotation invariant, then Gf(X.f,f t ) will be constant along geodesies for 
each X £ so(n). 

• The action of the scaling group of R given by c i— > e r /, with fundamental 
vector field £„(/) = a./. 

If the metric is scale invariant, then the scaling momentum G/(/, /t) will 
also be invariant along geodesies. 



5. The geodesic equation on Imm(M,R ra ) 



The G° metric is the simplest metric on Imm(Af , R n ). It is given by the following 
formula: 

G° f (h,k)= [ g(h,k)vol(g). 

This metric is well studied, see for example |15) . In this chapter we will study 
metrics of the form 



Gj(h,k)= / $(f)g(h(x),k(x))vo\(g)(x), 
Jm 

where $ is a Diff (M)-invariant function depending on the immersion / and possibly 
on x. These metrics are called almost local metrics. This definition includes as an 
important special case conformal versions of the G° metric, i.e. metrics of the form 



Jm 



where $ is again some Diff(M)-invariant function depending on the immersion / 
but not on x. Conformal metrics have been studied in [THl [T2]. We will consider 
functions $ depending on the volume and mean curvature, i.e. 

$ = $(Vol,Tr(i / )(a:)). 

For this class of functions, some work has already been done by [T7] for the special 
case of immersions of the unit circle in the plane. 
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5.1. The geodesic equation on Imm(A/, R"). We use the method of section 4.4 



to calculate the geodesic equation. So we need to compute the metric gradients. 
The calculation at the same time shows the existence of the gradients. Let m € 
T f lmm(M,R n ) with 



m = a.v^ + Tf.m 



T 



To shorten the notation, we will not always note the dependence on / in expressions 
as v> , TJ . . . 

D (f>m) Gj(h,k)= f (a 1 $)p (/ , m) Vol)5(/ l ,fc)vol(. 9 ) 



(d2$)(D if>m) Tr(L))g(h,k)vol(g) 
$.g(h,k)(D (f!m) vol(g)). 



To read off the i^-gradient of the metric, we write this expression as 



$.3 



M 



— (D u , m) Vol) + — (D {f . m) Tr(L)) + 



h, k I vol(g) 



Therefore, using the formulas from section [3] we can calculate the K gradient: 



Kf(m, h) 



— (£> (/>m) Vol) + — (£> (/>m) Tk(£)) + yol(g) 



ft* 



Tr(L).avol(g) 



+ ^-(-Aa + a Tr(i 2 ) + dTr(L)(m T )) 

+ div 9 (m T ) - Tr(L).al/i. 

To calculate the iJ-gradient, we treat the three summands of D^j m )Gj(h,k) sep- 
arately. The first summand is 



A/ 



(&*)(D (/im) Vol(a ; ))ff(/i(a ; ) ) vol(</)(x) 



xEM 



- Tr(L(y)).a(y) vo\(g)(y))g(h(x), k(x)) vo\(g)(x 



yeM 



g a(y).u(y),- Tr(L(y)) / (ft$)g(fc(aO, fc(x)) vol( ff )(x).^(y) vol( 3 )(y) 



xeM 



Gf(m, -^Tr(L) J (d^)g(h, k) vo\{g).v) . 



We will use the following formula, which is valid for any two functions a, b € 

C°°(M): 



A(a).bvol(g) 



M 



a.A(b)vo\(g). 



Thus the second summand is given by 

r (a 2 $)(£» (/>rn) Tr(i)) 3(/ l ,fc)vol( 5 ) 



/M 
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= / (<9 2 *)(- Aa + aTr(L 2 ) + dTr(L)(m T ))g(h,k)vo\(g) 
Jm 

= f -a.A((d 2 $)g(h,k))vol(g) + f a.(d 2 <S>)Tr(L 2 )g(h,k)vo\(g) 
Jm Jm 

+ [ (d 2 <f>)g( g ™d°(TY(L)),m T )g(h,k)vol(g) 
Jm 

= f -a.A((d 2 $)g(h,k))vol(g) + f a.(d 2 <S>)Tr(L 2 )g(h,k)vo\(g) 
Jm Jm 

+ I (9 2 $).9(T/.grad 9 (Tr( J L)),T/.m T )g(/i,fc)vol( ff ) 
Jm 

-G*(™,-^A((a 2 $).g(/,,fc)).^+G*(m,^(a 2 $)Tr(L 2 )g(/,,fc).^ 
+ G* (m, ^(d 2 $)g(h, k)Tf. grad^Tr(L))' 



In the calculation of the third term of the Hf(m, h) gradient, we will make use of 
the following formula, which is valid for e G°°(M) and X e X(M): 

0= f div(<jy.X).vo\(g) = [ £ . x vol(.g) 
Jm Jm 

= / (d o i^.x + i<p.x o d) vol(g) = / d{4>.i x vo\{g)) 
Jm Jm 

= / d(j) A ix vol(g) + / (j) A d(ix vol(g)) 
Jm Jm 

= / ( - ix(d<f>A vol(g)) +i x o d(p A vol(g)) + / 4>.C X vol(g) 
Jm Jm 

= I d<t>{X) vol(.g) + / <f>.div(X)vo\(g). 
Jm Jm 

Therefore we can calculate the third summand, which is given by 

/ *g(h, k)(D u , m) vol(g)) = [ <S>g(h, k)(div°(m T ) - Tr(L).a) vol(g) 
Jm Jm 

= - f d($g(h, k))(m T ) vol(g) + G*(m, -g(h, k) Tr(L).v) 
Jm 

= - J g(Tf. gmd 9 (<S>g(h, k)), ro) vol($) + G*(m, -g(h, k) Tr(L).v) 
= G* (m,-^Tf. gmd 9 ($g(h, k)) - g(h, k) Tr(L).//) 

Summing up all the terms the iJ-gradient is given by 

Hf{K k) = [ - | Tr(L) J (9i*)5(h, k) vol(g) - ±A((d 2 $)g(h, k)) 



M 

f 



^(8 2 <Z>) Tr(L 2 )g(h, k) - g(h, k) Tr(L) 



v 



+ 4> 



X -Tj. [(d 2 $)g(h, k) grad 9 (Tr(L)) - grad 9 (<i>^, k)) 
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Theorem. The geodesic equation for the almost local metrics G* on Imm(M, W 1 ) 
is then given by 



f t = h = a.vf + Tf.h T , 
1 



ht = 



i Tr(L) J \\h\\ 2 vol(.g) - ^A((9 2 $) \\h\\ 2 ) 



±(d 2 $)Tr(L 2 )\\hf -\\hfTr(L) 
1 

2$ 



Tf. ||/i|| 2 grad 9 (Tr(L)) - grad 9 ($ ||/i|| 



Tr(L).avol(f7) 



+ ^(-A fl + a Tr(L 2 ) + dTr(L){h T )) 
+ div 9 (h T ) -Tr(L).a]h. 



5.2. Momentum mappings. The metric G* is invariant under the action of the 
reparametrization group Diff(M) and under the Euclidean motion group M n xi 
SO(n). According to section [4. 5| the momentum mappings for these group actions 
are constant along any geodesic in Imm(M, 1"): 



VXeX(M): <f>(Vol(i),Tr(L f ))g(Tf.X,f t )vol(g) reparam. momenta 

J M 

$(Vol(f), Tr(L f ))g(fJ) vol(g) E T(T*M <S> M vol(M)) reparam. momentum 
/ <£(Vol(f), Ti{L* ))ft vol(g) linear momentum 

J M 

yXeso(n): / $(Vol(f), Tr(L f ))g(X.f, f t ) vol(g) angular momenta 
Jm 

angular momentum 



or I $(Vol(i),Tr(Z/))(/A/ 4 )vol(<?)€ A 2 R n = so(n) 

M 



6. The geodesic equation on Bi(M, R ra ) 

6.1. The horizontal bundle and the metric on the quotient space. Since 
vol(/*g) and Tr(L^) react equivariantly to the action of the group Diff(M), every 
G*-metric is Diff (M)-invariant. As described in Section 1.3 it induces a Riemann- 



ian metric on Bi (off the singularities) such that the projection 7r : 1mm — > Bi is a 
Riemannian submersion. 

By definition, a tangent vector h to / s Imm(M,K. n ) is horizontal if and only 
if it is G*-perpendicular to the Diff (M)-orbits. This is the case if and only if 
g(h(x),T x f.X x ) — at every point x <G M. Therefore the horizontal bundle at the 
point / equals the set of sections of the normal bundle (see Section 2.9 1 along / 
and thus the metric on the horizontal bundle is given by 



Gf(h hm ,k hm ) = Gj(a-v,b't/) 



$(Vol(£),Tr(L f )(x))a.bvol(g). 
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The following lemma shows that every path in Bi corresponds to exactly one 
horizontal path in Imm and therefore the calculation of the geodesic equation can 
be done on the horizontal bundle instead of on Bi. 

Lemma. For any smooth path f in Imm there exists a smooth path (p in Diff(M) 
with tp(t, ) — Idjv/ depending smoothly on f such that the path f(t,ip(t,x)) is 
horizonzal, i.e. dtf(t,ip(t,x)) lies in the horizontal bundle. 



The basic idea is to write the path if as the integral curve of a time dependent 
vector field. This method is called the Moser- Trick, (see [THl Section 2.5]) 



6.2. The geodesic equation on Bi(M, W 1 ). As described in section 1.3 geodesies 
in Bi correspond to horizontal geodesies in Imm. A horizontal geodesic / in Imm 
has ft = a.v^ with a g C°°(IR x M). The geodesic equation is given by 

ftt — CL t .v + a.v t — —H(a.v,a.v) — K(a.v,a.v), 

normal tang. 

see section [4. 4| This equation splits into a normal and a tangential part. From the 
conservation of the reparametrization momentum, see section 4.5 and the previous 
section, it follows that the tangential part of the geodesic equation is satisfied 
automatically. We will nevertheless check this by hand. From section [O] where 
we calculated the metric gradients on Imm, we get 

<9i$ 



Kf(a.v, av) 



Hf(a.is, av) 



Tr(L).avo\(g) 
- Aa + aTi(L 2 )) - Tr(L) .a] a.. 



Tf. 



(<9 2 $)a 2 grad 9 (Tr(L)) - grad 9 ($a 2 



Tr(L) / <9i$a 2 vol( ff ) - -A((<9 2 $)a 



M 



(9 2 $)Tr(L 2 )a 2 - a 2 Tr(i) 



From this we can easily read off the tangential part of the geodesic equation 
1 



a.v t 



2$ 



Tf. (a 2 $)a 2 grad 9 (Tr( J L)) -grad s ($a 2 ) 

To see that this equation is satisfied, we have to rewrite the right hand side using 
the following Leibnitz rule for the gradient, 

s(grad»(/i./ 2 ), X) = g{h grad 9 f 2 + f 2 grad 9 fi,X). 

This yields 

1 



a.v t 



2$ 
1 

2$ 
1 

2$ 



Tf. 
Tf. 



Tf. 



(<9 2 $)a 2 grad 9 (Tr(L)) - grad 9 ($a 2 
(<9 2 $)a 2 grad 9 (Tr(L)) - $. grad 9 (a 2 ) - a 2 . grad 9 ($) 
(<9 2 $)a 2 grad 9 (Tr(L)) - $.grad 9 (a 2 ) - (<9 2 $)a 2 grad 9 (Tr(L)) 



2$ 



$T/.grad 9 (a 2 ) = -Tf.a. grad 9 (a). 
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By the variational formula for v in section |3.11| this equation is satisfied automati- 
cally. The normal part is given by 

a t — g{ -^B(a.v, a - v ) ~ K(a.v, a.v), v) 

= i[^a 2 Tr(L/) -±Tz{lJ) J ™Kf*g) - ^,A((d 2 $).a 2 ) 

+ (Si$)a / Tr(L- f ).avol(/*g) - Tr((Z/)V + (d 2 $)aAa\ . 

We can rewrite this equation by expanding Laplacians of products as follows: 

A (0102) = (Aoi)oa - 2Tr 9 (cZai da 2 ) + ai(Aa 2 ). 
Theorem. The geodesic equation of the almost local metric G* on Bi reads as 
ft = a.v f , 

at = ^[^« 2 Tr(L0 - ^ Tr(^) ^(^^a 2 vol(/* .g) - i a 2 A(3 2 $) 
+ 2aTr 9 (rf(<9 2 $) <g) da) + (<9 2 $) Tr s (da ® da) 
+ (^$)a / Tr(i/).avol(rg)-i(9 2 $)Tr((L0 2 )a 2 ]. 

For the case of curves immersed in K 2 , this formula specializes to the formula 
given in [TTJ section 3.4]. (When verifying this, remember that A = — D 2 in the 
notation of [17].) 

7. Sectional curvature on shape space 

To compute the sectional curvature we will use the following formula, which is 
valid in a chart: 

#0(01,02,01,03) = Gq (i2o(ai, 02)01, 02) = 

^d 2 G* (ai,Oi)(a 2 ,02) - d 2 G* (ai,a 2 )(ai,a 2 ) + ^d 2 G* (a 2 , a 2 )(ai, ai) 

+ G* (r (ai,ai),r (a 2 ,a 2 )) - G* (T (ai , a 2 ) , T (ai , a 2 ) ) . 

Sectional curvature is given by 

#0(01,02,02,01) = -^0(01,02,01,02). 

Therefore we have to calculate the metric in a chart and calculate its second deriv- 
ative. 

7.1. The almost local metric G* in a chart. In the following section we will 
follow the method of [T5]. First we will construct a local chart for Bi. Let /o : 
M — > R™ be a fixed immersion, which will be the center of our chart. Consider the 
mapping 

ip = tp fo : G°°(M, (-e, e)) Imm(M,K") 
^{a){x) = e^% (x) (a(x).vf°{x)) = f (x) + a{x).v f °(x), 
where e is so small that ?p(a) is an immersion for each a. 
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Denote by w the projection from Imm(M,R") to Bj(M,R n ). The inverse on its 
image of tt o V> : C°° (M, (-e, e)) ^ B, (M, R") is then a smooth chart on (M, R"). 
We want to calculate the induced metric in this chart, i.e. 

((7roV0*G*)a(&l,&2) 

for any a € C°°(M, (-e,e)) and 61, 6 2 G C°°(M). We shall fix the function a and 
work with the ray of points t.a in this chart. Everything will revolve around the 
map: 

f(t, x) = i){t.a){x) = f (x) + t.a(x).v /o (x) 

We shall use a fixed chart (u,U) on M with di — ^j. Then in this chart, the 
pullback metric is given by 

n— 1 n— 1 

g\u = ^2 9t du l <g> du 1 = ^2 9{dif, djfidu 1 <g> dw 3 , 

the volume density by 

vol(g) = y/det(g {dj^dj])^ 1 A • • • A du"" 1 !, 
the second fundamental form by 



,H - 8{d it dj) = gWa.Tf.dj, v) = g[^l v h), 

and the mean curvature by Tr(L) = J2i j g K ' s ij- To calculate the metric G* in this 
chart we have to understand how 

rt. V.6i =b 1 (x).v fo (x) 

splits into a tangential and horizontal part with respect to the immersion f(t, ). 
The tangential part locally has the form 

n-l 

T/.(T (t . a) ^.(6i)) T = ^c^ J /(t,x), 
where the coefficients c l are given by 

n-l 

= ^2g ii g(bi(xyo{ x ),d j f{t,x)) . 

3 = 1 

Thus the horizontal part is 

n-l 

{Tt.^.b^.u = (Tt.^.h) - T/.(T (t . a) ^.(6 1 )) T = h(x)^(x) - £ c%f(t,x). 

»=i 

Lemma. TTie expression o/G* m i/ie c/iart (ttot/;) -1 is; 

((TT O V/„)*G*)( t .a)(6l, 62) = G? w(t . o)) (T t . (7T O ^).6 1 ,T t . (7T o V>)6 2 ) 

= G| (4 . a) ((T t . a ^.6 1 ) ± .^ (^.62)^.^) 

$5 ((Tt^V-fci)" 1 .!/, (Tt.aV.63)) vol(. 9 ) 



/ 

J A 



M 

n-l 



= \ b 1 .b 2 -J2^9(dif(t,x)Mx)-v ,0 {x))J vol(g) 
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7.2. Second derivative of the G*-metric in the chart. We will calculate 

a t 2 |o((^o^ )*G*) (t . a) (6 l7 6 2 ). 
We will use the following arguments repeatedly: 

d^f = djd t \ f = = (d ja y° + a(djvf°), 

tang. 

^(x^x),^/^)) |t=o = 0, 
and consequently Cj| t =o = 0. 

n— 1 n— 1 

ftloc, = X>|o(0 ij > + X> tf ft|off (6i*>,^/) 

3=1 3 = 1 

n—1 n—1 7i — l 

= E ^ (6i^ /o , a t | a,/) - Yl 9 ij S (M A , d,(a.z>)) = 
i=i i=i j=i 

Therefore 

n-l 

(6i.6 2 -Y, ci 9(dif, b 2 .y fa )) | t=0 = 61.62 



»=i 

n-l 



9 t |o(6i.62-Ec^(ft/,62^ /o )) = 

i=l 

n—1 n—1 

= -J^lo^O- E°- 9 *lo^/,62.^ /o ) = 

i=l i=l 
n-l 

9 t 2 |o(6i.62-Ec^/,6 2 .^ )) = 

n—1 n—1 n—1 

= -E(a 2 | oC l ).0-2E(a t |oc l )9 t |o5(aj,62.^°)-E0.a t 2 |o.9(^/,6 2 .^) 

-i— 1 i — l i—1 

n—1 n—1 

= -2Y / (dt\^)g(d l (a.^),b 2 .^) = -2^(9 t | oC ')(9 lS )6 2 

i=l i=l 
n— 1 n— 1 

= -26i6 2 Yl = - 2& i fo 2 ||da||J-i • 
»=i j=i 

The derivatives of $ are 

a t | ($o (Vol,Tr(L))) = (9 1 $).(a t | Vol) + (9 2 $).(9 t | Tr(L)) 

9 2 |o($o (Vol,Tr(L))) = (9 1 a 1 $).(a t | Vol) 2 + (a 2 a 2 $).(a t | Tr( J L)) 2 

+ 2(9 1 a 2 $).(9 t | Vol).(5 t | Tr(L)) + (Si$)(d t 2 | Vol) + (d 2 $)(d 2 | Tr(L)). 

Lemma. TTie second derivative of the G* -metric in the chart (ir o is groen 

5 2 |o((7ro^ /o )*G*) (ta) (6i,6 2 ) = (d 2 ((7ro^ /o )*G*)(0)(a,a))(6i,6 2 ) 
(1) /■ 

= / ... 61.62 VOl(fif) 
JM 
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over the following expression 

.-. = - 2 HifaH^ ) + (ft*).((#|o Vol) + 2(ft| Vol)^) 

+ (5 2 $).((9 t 2 | Tr(L)) + 2(9,| Tr(i))^fc 01 ) + (d^).^ Vol) 2 
+ 2(a 1 9 2 $).(a t | Vol)(9 t | Tr(L)) + (a 2 a 2 <i>).(a f | Tr(L)) 2 . 

7.3. Sectional curvature on shape space. To understand the structure in the 
formulas for the sectional curvature tensor, we will use some facts from linear 
algebra. 

7.3.1. Lemma. Let V = C°°(M), and let P and Q be bilinear and symmetric maps 

V x V -> V. Then 

ffl(P, Q)(ai A 02, h A 6a) := ^(P{a 1 ,b 1 )Q{a 2 ,b 2 ) - P(a 1 ,b 2 )Q{a 2 , h) 

+ P(a 2 , b 2 )Q{a 1 ,b 1 ) - P{a 2 ,b 1 )Q{a 1 ,b 2 )) 
defines a symmetric, bilinear map (V A V) ® (V A V) — > V . 

Also EB(P, Q) = EH(<2, -P). The symbol EE stands for the Young tableau encoding 
the symmetries, see [7J. We have 

ffl(P,Q)(oi A a 2 ,a 1 A a 2 ) 

= -P(a 1 ,a 1 )Q(a 2 ,a 2 ) - P(ai, a 2 )Q(a 2 , m) + -P(a 2 , a 2 )<3(ai, ai). 

P is called positive semidefinite if for all x € M and a £ C°°(M), P(a, a)(x) > 0. 
P is called negative semidefmite if — P is positive semidefmite. We will write P > 
0, P < 0, P > if P is positive semidefmite, negative semidefmite or indefinite. 

7.3.2. Lemma. If P and Q are positive semidefinite bilinear and symmetric maps 

V x V —¥ V , then also EB(P, Q) is a positive semidefinite symmetric, bilinear map. 

Proof. To shorten notation, we will write for instance P\ 2 instead of P{a\, a 2 ). The 
Cauchy inequality applied to P and Q gives us 

P12Q12 < yj P11P22Q11Q22, 

and therefore we have 

EE!(P, Q)(oi A a 2 , ax A a 2 ) = -PuQu - P12Q12 + \P22Q22 

> -^PwQii — yj P11P22Q11Q22 + 2-F22Q11 

= ^(^11^22- VP22Qll) 2 >0. □ 

Let A, ^ : V — > V. Then the map A®/i:V®V— ^Vis given by 

(A <x> /i)(a ® 6) = X(a).fi(b), 
where the multiplication is in V = C°° (M) . Denote by A V /1 the symmetrization 
of the tensor product given by 

AV//= §(A® /x + jt® A). 
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We will make use of the following simplifications: 

7.3.3. Lemma. Let A, /?, fi, v, : V V . Then the bilinear symetric map 

ffl(A V^,/iV v) 

satisfies the following properties: 

(51) EH (A V (j,, A V v)(ax A a 2 , a x A a 2 ) = — ~(A ® n){ct\ A a 2 ).(A <g v)(ai A a 2 ), 

(52) ffl(AV/i,A®A) = 0, 

(53) EB(A ® A,/x V i/)(ai A a 2 ,ai Aa 2 ) = ^(A® A 4 ) (ax A a 2 ).(A ® f)(ffli A a 2 ). 



Proof. For the proof of simplification (SI ) we calculate 



3(A V /i, A V f)(cii A a 2 , ai A a 2 ) 
f 



-(A ® |U ® A ® v) 





$ a 2 <i 


§ a 2 + 


a 2 $ 


3 a 2 <5 


5 ai G 


§ ax 


a 2 <5 


5 ai 5 


S a 2 - 


f 


5 a 2 <5 


5 a 2 $ 


5 ax 


ai <i 


5 Oi 5 


5 a 2 - 


f 

-a 2 d 


5 ai $ 


5 a 2 $ 


5 ax 



Using the symmetries of the quasilinear mapping A ® /i <E> A <g /i, we can swap the 
first and third position in the tensor product of the two summands in the first line. 
Then the expression inside the square brackets equals — \{a\ A a 2 ) <g (ax A a 2 ). 



Since A ® A vanishes when applied to elements of V A V, simplification (S2 ) is a 



direct consequence of (SI I 



For the proof of simplification (S3 1 we calculate: 



3(A ® A, ^ V v)(ai A a 2 , a>i A a 2 ) 
1 



-(A® A® /u® i/) 



ax ® ax ® a 2 <g a 2 + a 2 <g a 2 ® ax ® ax 
ax (g a 2 (g ax (g a 2 — ax (g a 2 (g a 2 (g ax 



Using symmetries as above, we can replace third summand a x ® a 2 ® ax (g a 2 by 
a 2 ® ax (g a 2 ® ax, because the first two tensor components oi\®\®^®v are 
equal. Then, swapping the second and third position in all tensor products, we get 

EH (A ® A, fi ® v)(a\ A a 2 , ax A a 2 ) 

= - (A ® /i ® A ® v) ai <E> a 2 <E> a x <E> a 2 + a 2 <gi a x <E> a 2 <E> ai 
—a 2 <g ax <g a x <g a 2 — ax <g a 2 <g a 2 <g ax 
The expression inside the square brackets equals (ax A a 2 ) ® (ax A a 2 ). □ 



For orthonormal ax, a 2 sectional curvature is the negative of the curvature tensor 
i?o(ax, a 2 , a 1; a 2 ). We will use the following formula for the curvature tensor, which 
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is valid in a chart: 

Ro(ai,a 2 ,ai,a 2 ) = G® (R (ai, a 2 )ai, a 2 ) 
1 



(1 



-d G (ai, Oi)(oa, a 2 ) - d G (oi, a 2 )(ai, a 2 ) + -d G (a 2 , a 2 )(ai, ai) 
+ (r (ai,ai),r (a 2 , a 2 )) - G* (r (ai, a 2 ), r (ai, a 2 )). 



Looking at Formula ([I]) from section 7.2 we can express the second derivative of 
the metric G* in the chart as 

d 2 (^o^ /o )*G*(0)(a 1 ,a 2 ))(6 1 ,6 2 ) 

fs.PrKaa) + (a 1 $)P 2 (a 1 ,a 2 ) + («9 2 $)P 3 (ai, a 2 ) + (5 1 9 1 $)P 4 (a 1 , a 2 ) 

A/ V 

+ (3id 2 $)P 5 (ai, 02) + (9 2 a 2 $)P 6 (ai, a 2 ))p(6i, & 2 ) vol(g), 

where P{b\,b 2 ) — b±.b 2 , so P = id® id, and where the Pj are obtained by sym- 
metrizing the terms in Formula ([!]) from section 7.2 



For the rest of this section, we do not note the pullback via the chart anymore, 
writing G* instead of ((vr o V/ )*G*) (0), for example. To further shorten our 
notation, we write L instead of and g instead of go. The following terms are 
calculated using the variational formulas from section [3] 



Pi (a, a) 



P 4 (a, a) 
P 5 (a, a) 



cfjovol 
vol 



P 2 (a,a) = (#|oVol) + 2(&| Vol) 



2||da||^ 1 =a 2 (Tr(L) 2 -Tr(L 2 ))-||da||^_ 1 
3J vol 



vol 



a 2 (Tr(L) 2 -Tr(L 2 )) + / \\da\\^ vol( 5 ) 

A/ JM 



+ 2Tr(L).a / Tr(i).a vol(g) 
Ja/ 

P 3 (a,a) = (3 2 | Tr(L)) + 2 (ft| T^)^p 

= 2a 2 Tr(L 3 ) + 4a Tr (L.0 -1 .V 2 a) + 2Tr(g _1 (da (g) da)L) 

- \\da\\ 2 g ^ Tr(L) + 2aTr s (V gradaS ) 

+ 2(- Aa + aTr(P 2 ))(-Tr(P).a) 
= 2a 2 Tr(L 3 ) +4aTr( J L..g" 1 .V 2 a) + 2 Tr^" 1 (da ® da)L) 

-||da|| 2 _ 1 Tr(L) + 2aTr 9 (V 

grad a$) 

+ 2 Tr(L)aAa - 2 Tr(i) Tr(L 2 ).a 2 
(ftloVol) 2 = f / Tr(L).avol(.g)" 2 

M 



2(9 t |oVol)(9 t | Tr(£)) =2 / - Tr(L).avol( 5 )( - Aa + aTr(L 2 )) 

JAf 

2Aa / Tr(i).avol(g) - 2Tr(L 2 )a / Tr(L).a vol(p) 

JM JM 
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P 6 (a,a) = (d t \ Tr(L)) 2 = ( - Aa + aTr(Z 2 )) 2 
= (Aa) 2 - 2aAa Tr(L 2 ) + a 2 Tr(L 2 ) 2 
Then the first part of the curvature tensor is given by 

^d 2 G* (ai, ai)(a 2 , a%) - d 2 G$ (ai, a 2 )(ai, a 2 ) + ^cPG* (a 2 , a 2 )(ai, oi) 

= / ($.ffl(P 1 ,P) + (9 1 $)ffl(P 2 ,P) + (9 2 $)ffl(F 3 ,P) 

J AI 

+ (dxdx®) ffl(P 4 , P) + (9i9 2 $) ffl(P 5 , P) + (d 2 d 2 $) ffl(P 6 , P)) vol(. 9 ) 
• (ai A a 2 , ai A a 2 ). 

Note that P is positive definite, so that EB(Pi, P) is positive semidcfinite if Pi is 
positive semidcfinite. We can always assume that $ is positive because otherwise 
G* would not be a Ricmannian metric. 



PP 



with 



P =P 1 1 +P 2 , 

P^ = (Tr(P) 2 - Tr(L 2 ))id®id 
P 2 = _Tr 9 (d(g)d) 



Applying simplification ([S3]) to ffl(P 1 1 ,P) and ffl(P 1 2 ,P), we get 



i{P*,P) = -(Tr(L) 2 - Tr(L 2 ))(id®id) 2 = 



on (V A V) ® (7 A V) and 



3(P 2 ,P) = --Tr2((id^) 2 ), 

K-ff 5 - p )( a i A fl 2,«i A a 2 ) = -~ ||ai«fa 2 - » 2 d»i|| 2 -i < 0. 



Therefore we have 



P 2 P 



with 



/ $.ffl(Pi,P)(ai Aa 2 ,cti Aa 2 )vo%) < 0. 



P 2 = Po 1 + P| + P| 



P 2 : = / (id®id)(Tr(L) 2 - Tr(L 2 )) vol( 5 ) 
P| = 2 Tr(X) (id V / Tr(L)idvol(#)) 
Tr 9 (d<g)d)vol(g) 



P? = 



A/ 
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P 2 is indefinite. Applying Simplification ( S2 1 we get EH (Pi, P) — 0. P 2 an d therefore 
also EH(P|,P) is positive semidefinite. Therefore 



/ (ai$)ffl(P 2 x ,P)(ai Aa 2 ,a 1 Aa 2 )vol(<?) < 0, 

/ (<9i$) ffl(P|,P)(ai A a 2 , ai A a 2 ) vol(fif) > 
Jm 



PzP 



Pa = p} + P? + P? 



3 ) 



P3 1 = 2idV Tr(L 3 )id+2Tr(L<r 1 V 2 (id)) + Tr s (V gradld s 



with 



+ Tr(i)A(id) - Tr(i) Tr(L) 2 id 
P| = 2Tr(g" 1 (d<g>d)L) 
P 3 3 = -Tr s (d(g)d)Tr(i) 

Applying Simplification ( |S2[ ) we get that EH(Pg,P) vanishes. Furthermore, 

ffl(P 3 2 ,P)(ai A a 2 , ai A a 2 ) = a\ Tr(g~ 1 (da 2 <g> da 2 ).L) 
— 2aia 2 Tr^ -1 ^^ (g) da 2 ).L) 
+ a 2 Tr(cT :L ((iai <8>dai).P) 
= g 2 (( a i^ a 2 — a 2 c£ai) (g) (airfa 2 — a 2 <iai), s)) > 



3(P 3 3 , P)(ai A a 2 , 01 A a 2 ) = -- ||ai<fa 2 - flfcdaiHo-i Tr(£) < 



P*P 



Pa= { Tr(i)idvol(5)(g) \ Tr(L) id vol(g) 
Jm Jm 



Applying Simplification ( S3 ) we get 



EH(P 4 ,P) = ~(i,l 



Tr(L) idvol(s) 



M 



Therefore, if <9i<9i$ > 

/ 9 1 9i$ffl(P 4 ,P)(ai Aa 2 ,ai Aa 2 )vol( 5 ) > 0, 
Jm 



PnP 



M 

P 5 = P 5 1 +P| 



with 



Pi = 2 (A V / Tr(i) id vol(g) 
^ Jm 

Pf = -2 Tr(L 2 )( id V / Tr(i)a 2 vol(.g) 



Applying Simplification (S3) we get that ffl(P 5 1 ,P) is the indefinite form given by 

ffl(P 5 \P) = (id(g>A) ® (id® f Tr(L)idvol( ff )) 

Jm 
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Simplification ([S2} gives ffl(P 5 2 ,P) = 0. Therefore 

/ (d 1 d 2 ^)m{P^P){a 1 Aa 2 ,a 1 Aa 2 )vo\(g) < 0. 



P fi P 



P 6 = P, 1 + Pi 



with 



P 6 = A (g) A 



1 

P 6 2 = Tr(L 2 ) 2 id®id 
P| = -2Tr(L 2 )idVA 



Applying Simplification (S2) we get that ffl(P|,P) and ffl(Pg,P) vanish. Simplifi- 
cation (S3 1 gives 



3(P 6 \P)= 2 (id®A) 2 



We get 



/ (<9 2 a 2 $)ffl(P 6 ,P)( ai Aa 2 , ai Aa 2 )vol( 5 ) > 

J M 



if a 2 9 2 $ > o. 



Now we come to the second part of the curvature tensor R${ai, a 2l a\, a 2 ), which 
is given by 

G (r (ai, ai), r (a 2 , a 2 )) - G (T (ai, 02), T (ai, a 2 )). 
From the geodesic equation calculated in section [6j which is given by 

at = T (a,a) = -j- [i$a 2 Tr(L) - ^ Tr(L) J^(di$)a 2 vo\(g) - \a 2 A(d 2 <S>) 

+ 2aTr 5 (d(<9 2 $) <g) da) + (9 2 $) Tr 9 (da ® da) 

+ (d x $)a / Tr(L).avolQ?) - i(9 2 $) Tr(L 2 )a 2 ] , 

we can extract the Christoffel symbol by symmetrization and get 

5 



r (ai,a 2 ) = ^r^Q^a^aa), 



where Qi, . . . , Q5 are the symmetrizations of the summands in the geodesic equa- 
tion. Qi are given by 

Qi = \{® Tr(L) - A(9 2 $) - (d 2 $) Tr(L 2 )) id ® id, 
Q 2 = -\ti{L) f (c\$)id®idvol( ff ), 

* J M 

Q 3 = 2idVTr s (d(<9 2 $)<g>d) 

Qi = (9i$)idV / Tr(i)idvo%), 
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Q 5 = (d 2 $)Tr 9 (d(g>d). 



Then 



G (r (ai, ai), r (a 2 , a 2 )) - G (T (ai, a 2 ), T (ai, a 2 )) 
= I *y2 m (QhQi)( a l A «2,«i Aa 2 )vol(5) 

Jm * ; 

+ / ty2&(Qi,Qj)(ai Aa 2 ,ai Aa 2 )vol(g). 



l<3 



The contribution of the following terms to i?o( a i, a 2, ai, a 2) is J M 5 ■ ■ ■ vol (5) 
over the terms listed. 



Q1Q1 ffl(Qi,Qi) = 



according to Simplification ( S2 1 



Q2Q2 ffl(Q 2 , Q 2 )(ai A a 2 ,ai A a 2 ) = 



Tr(L) 5 



(5 1 $)a?vol(«?). / (5 1 $)a|vol( 9 )-( / (di$Ka 2 vol(<?))' 
which is positive by the Cauchy-Schwarz inequality, assuming that di<£> > 0. 



Q3Q3I ffl(Q3, Qz)(ai Aa 2 ,fli A a 2 ) = 

= -((id(8)Tr 9 (d(a 2 $)®d))( ai Aa 2 ) 
= — .g _1 (d(9 2 <i>),ai(ia 2 — a 2 dai) 2 < 



according to Simplification ( |S1[ ). 

Q4Q4I ffl(Q 4 ,Q 4 ) = -^(9i$) 2 (id® / Tr(L)idvol( 5 )) 2 < 
4 Jm 



according to Simplification (SI I. 



Q5Q5 ffl(Q 5 ,Q 5 ) = (d 2 $) 2 (||d ai || 2 x \\da 2 \\ 2 ! -^(doijdoa) 2 ) 



(a 2 $) 2 ||d ai Ada 2 ||gO > 



by the Cauchy-Schwarz inequality. 



The contribution of the following terms to i?o(ai, 12,01,12) is J M 5 • ■ ■ vol(g) over 
the terms listed. 



QiQ 2 | ffl(Qi,Q 2 ) - -|($Tr(L) 2 - Tr(L)A(d 2 $) -Tr(L)Tr(L 2 )(d 2 $)). 



(id<8>id, / (ai$)id(8)idvol(3)V 
^ Jm ' 
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where the second factor is > assuming that 9i$ > 0. 
\Ch\ ffl(Qx,Q 3 )=0 



according to Simplification (S2|. 
QiQA ffl(Qi,Q 4 ) = 



according to Simplification (IS2 



1 



QiQi 



H(Qi, Q 5 ) = ~ ($ Tr(L)(0 2 $) - (0 2 $)A(0 2 $) - Tr(L 2 )(0 2 $) 2 ) 



|aida2 — a 2 dai|| s _! 



Q 2 Qs 



ffl(Q 2 ,Q 3 ) £ 



Q2Q4 



m(Q 2 ,Q 4 ) = --(d 1 ^)Tr(L)- 



M 

This form is indefinite, but we have 
2 



(<9!<I>)id«)idvol(3),idV / Tr(L) id volte)) 



A/ 



A/ 



3(Q 2 ,Q 4 ) vol(. 9 ) 



with the positive semidefinite form 

(9i$)id(g)idvol(g), 



A/ 



and the form 



1 



"&(£)- (01$) id vol(3) V / Tr(i) id volte), 



A/ 



A/ 



which is positive semidefinite if ^j- is a non-negative constant. 



Q2Q1 



S(Q 2 ,Q 5 ) = --(d 2 ^)Tr(L)- 



A/ 



(0 1 $)id®idvol( 5 ),Tr ff (d(g)d)) < 0, 



because of the factor (0 2 $)Tr(L). But the factor 
(0!$) id® id vol( ff ), Tr ff (d (g) cf)) 



M 



is positive definite. 



Q3Q4 ffl(Q 3 ,Q 4 )^0 

ffl(<?3, Q5)(ai A a 2 , ai A a 2 ) = 

= (0 2 $)(ai5- 1 te(0 2 $),dai) HdasllJ-! - (a^ 1 {d{d 2 <S>) , da 2 ) 
a 2 g _1 (d(0 2 $), dai))g _1 (dai, da 2 ) + a 2 g' 1 (d(d 2 ^),da 2 ) \\dai 



+ 
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QiQl 



= {82^)9° (d(d 2 $) <8> {aida 2 ~ a 2 da 1 ),dai Ada 2 ) > 

ffl(Q 4 ,Q 5 ) £ 



We are now able to compile a list of all negative, positive and indefinite terms of 
Ro(ai, a%, <zi, a%). Remember that negative terms of i?o( a l> a 2> a li a %) make a pos- 
itive contribution to sectional curvature. Positive sectional curvature is connected 
to the vanishing of geodesic distance because the space wraps up on itself in tighter 
and tighter ways. 



P 4 P 


p^p 


Q2Q2 


Q5Q5 




Pip 


Q3Q3 


Q4Q4 


Q1Q2 



Q2Q2 Q5Q5 are positive, assuming <9i$, didi$, 8282$ > 0. 



Q2Q4 is negative assuming that is a non- negative constant, and indefinite 



otherwise. 



Q2Q5 is negative assuming that Tr(L)(9 2 ( I ) ) is positive, and indefinite otherwise. 



P2P 


P3P 


P 5 P 


Q1Q5 


Q2Q3 


Q3Q4 


Q3Q5 


Q4Q5 



8. Geodesic distance on Bi(M,R n ) 

We will state some conditions on $ ensuring that the almost local metric G* 
induces non- vanishing geodesic distance on Bi. The proofs are based on a compar- 
ison between the G*-length of a path and its area swept out. In the last part we 
will use the vector space structure of R" to define a Frechet metric on shape space 



Bi(M, W 1 ). In 8.7 it is shown how this metric is related to a L°° Finsler metric, 



and in |8.8| the Frechet metric is compared to almost local metrics. 
The main results are in section 18.61 and in section l8T 



8.1. Geodesic distance on Bi. Geodesic distance on Bi is given by 

distJ(/ ,/i)=infLg*(/), 

where the infimum is taken over all / : [0, 1] — > Bi with /(0) = /o and /(l) = f\. 
is the length of paths in Bi given by 



L%{f) 



G*(f t ,ft)dt for / : [0, 1] -> B t 



Letting ir : Imm — > Bi denote the projection, we have 
L ( T(irof)=L£ m (f) 



G?(/t, ft)dt for horizontal / : [0, 1] -> Imm. 



By non-vanishing geodesic distance on Bi we mean that dist B . separates points. 
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8.2. Area swept out. For / : [0, 1] — > Imm we have 

(area swept out by /) = / vol(/(-, ■)*<?). 

J[0,l]xM 

If / is horizontal, then this integral can be rewritten as 

(area swept out by /) = f f ||/ t || vol(/(t, -f^dt =: [ [ \\ft\\ vol(g)dt, 

JO JM JQ JAI 

8.3. Lemma (First area swept out bound). For an almost local metric G* satisfying 

$(Vol, Tr(I/)) > Ci for d > 0. 
and a horizontal path f : [0, 1] — > Imm, we have the area swept out bound 



\/~C\ (area swept out by f) < max -y/Vol (f(t)) -^iL m (/) 

The proof is an adaptation of the one given in jT5J section 3.4] for the G A -metric. 
Proof. 

L?lM) = J yjGj(f t ,ft)dt 

(f <f>\\f t \\\ol(g)) l2 dt>^[([ \\f t \\ 2 voUg)) l2 dt 
o k Jm ' Jo k Jm 

i 



> VCI / f / vol(. g )) 5 / l.||/ t ||vo%)di 
Jo k Jm ' Jm 

>V^min(/ vo\(g)Y 2 f f 1. ||/ t || vol(g)dt 
' k Jm ' Jo Jm 
1 pi 



= ^(max/ vol(g)) '/ / 1. ||/ t || vol(g)dt □ 
v ' Jm ' Jo Jm 

8.4. Lemma (Lipschitz continuity of VVol) . For an almost local metric G*, the 
condition 

$(Vol,Tr(L)) > G 2 Tr(L) 2 
implies the Lipschitz continuity of the map 

VVo\ : [B u distal,) K> 
fry £/ie inequality holding for f\ and fi in Bi : 

The proof is an adaptation of the one given in [IS] section 3.3] for the G A -metric. 

Proof. Let / : [0, 1] — > Imm be a horizontal path, and let f t — a.v^ denote its 
derivative. Using the formula from section 3.7 for the variation of the volume we 

get 

9 t Vol(/) = - / Tr(L)avo%)< / Tr(L)avol( 5 ) 
Jm Jm 
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<(/ l 2 vol(g)) 5 (/ Tr(L)Vvol( 5 )) 5 

< v^ol(7)( J ^a 2 vol(g))' < ly/vdltfjy/GjVtJt). 



Thus 



2v/Vol(7) " 2 ^ 
By integration we get 

v/voiC/iT- v / Voi(^)= / d t ^JvdU)dt 



(i 



1 i i 



, 2VC 2 V J """" 2^G 2 
Now take the infimum over all horizontal paths / connecting /i and / 2 . □ 

8.5. Lemma (Second area swept out bound). For an almost local metric G* sat- 
isfying 

$(Vol, Tr(L)) > G Vol wii/i C > 
and a horizontal path f : [0, 1] — > Imm, we gei the area swept out bound 

VC (area swept out by f) < ii^ m (/), 

The proof is adapted from proofs for the case of planar curves that can be found 
in [TTJ section 3.7], jTHl Lemma 3.2], proposition 1] and theorem 7.5]. 

Proof. 



L?L(f)= J sjGj(f t Jt)dt= J (y^$||/ f || 2 vol(.g)) 5 dt 

>VC [ VVoUJ)( [ ||/J 2 vol(g)) 5 di> VC f f l.||/ t ||vol( ff )di 
Jo k Jm ' Jo Jm 

= VC [ vol(/(-, -)*g)dt = \fC (area swept out by /). □ 

J\0,l]xM 



8.6. Non-vanishing geodesic distance. Using the estimates proven above, we 
get the following result: 

Theorem. At least on B e , the almost local metric G* induces non-vanishing geo- 
desic distance if one of the two following conditions holds: 

(1) $(Vol, Tr(£)) > Gi + G 2 Tr(L) 2 for Gi, G 2 > 0. 

(2) $(Vol, Tr(£)) > G Vol for G > 0. 
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8.7. Frechet distance and Finsler metric. The Frechet distance on shape space 
Bi(M, R n ) is defined as 

di S t%(F ,F 1 )= inf H/o-ML-, 

J0,J1 

where the infimum is taken over all /o,/i with 7r(/o) = i<o,7r(/i) = F\. As before, 
7r denotes the projection 7r : Imm — > . Fixing /o and /i , one has 

dist^(7r(/ ),7r(/i)) = inf ||/ o ip - /i|| ioo , 

where the infimum is taken over all cp G Diff(M). The Frechet distance is related 
to the Finsler metric 

G°° : rimm(M,M n ) — > M, ft, i-» ||/i- L || £oo . 

Lemma. TTie patft length distance induced by the Finsler metric G°° provides an 
upper bound for the Frechet distance: 



distS (F , F x ) < distgtc (F , fi) = inf / 1 1 /, 



f II qoo 



where the infimum is taken over all paths 

/: [0,1] ->Imm(M,R n ) witA tt(/(0)) = F , 4/(1))=^. 

Proof. Since any path / can be reparametrized such that ft is normal to /, one has 

f 1 /•! 



J Jo ' J Jo 



where the infimum is taken over the same class of paths / as described above. 
Therefore 



dist£'(F , Fi) = inf ||/(1) - /(0)|| z . = inf 



f ftdt <inf / ||/ t || £ao dt 

JO L~ J JO 



inf / ll/^H^ dt = dist^i o (F ,F 1 ), a 



It is claimed in [12l theorem 13] that doo = dist^no. However, the proof given 
there only works on the vector space C°°(M,M n ) and not on Bi(M, M"). The reason 
is that convex combinations of immersions are used in the proof, but that the space 
of immersions is not convex. 

8.8. Theorem (Almost local versus Frechet distance on shape space). On the shape 
space Bi{M, W 1 ) the G* distance can not be bounded from below by the Frechet 
distance if one of the following conditions holds: 

(1) $ < C x + C 2 Tr(i) 2fe for Cx,C 2 > and k < (dim(M) + 2)/2, 

(2) $<CVol ft for C>0, 

(3) $ < Ce y ° ! for OQ, 

(4) $ < Cx + C 2 de\(L) 21 for Cx,C 2 > and I < ~ ' 



2 dim(M) 



Indeed, then the identity map 

Id : (Bi(M,R n ),d G *) -> (A(M,R n ) )C U 
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is not continuous. 



Proof. Let /o be a fixed immersion of M into R™, and let f\ be a translation of 
fo by a vector h of length I. We will show that the IP-distance between 7r(/o) 
and 7r(/i) is bounded by a constant 2L that does not depend on I. It follows that 
the 7? p -distance can not be bounded from below by the Frechet distance, and this 
proves the claim. 

For small r , we calculate the G*-length of the following path of immersions: 
First scale fo by a factor ro, then translate it by h, and then scale it again until 
it has reached f\. The following calculation shows that under one of the above 
assumption the immersion fo can be scaled down to zero in finite G*-path length 
L. 



L lm m ( r -fo) = 




$(Vo\(r.fo),Tr(Lrfo),det(Lrf«))g(r t .f 0: r t .f )vol((r.f )*g)dt 

M 

f 1 J f rl<S>( r ™ Vol(/ ), - Tv(Lfo), -L det(Lf°))g(f , / )r™ vol (f*g)dt 
Jo v Jm r r 

$(f"Vol(/ ), -Tr(L/o), ±det(Lf°))g(f ,fo)r m vol(f5g)dr 

M r r 

The last integral converges for all of the above assumptions. Scaling down to r > 
needs even less effort. So we see that the length of the shrinking and growing part 
of the path is bounded by 2L. 

The energy needed for a pure translation of the scaled immersion by distance £ 
is given by (/ = t.h, with g(h, h) = I 2 ): 

L?L{f)= I I $.g(ftJt)vol(g)dt 



Jo Jm 

[ [ $.t 2 £ 2 vol{g)dt = £ 2 [ $vol(c?) 
Jo Jm Jm 
' 0( Y(m-2fc)) ; if $ satisfics g 

0(r m ( k + 1 )) 1 if $ satisfies ^ 

0(e r m .r m ), if $ satisfies ^ 

This length tends to zero as r tends to zero. Therefore 

distg* (tt(/o), tt(/i)) < dist^ m (/ , /i) < 2L. □ 



9. The set of concentric spheres 



For an almost local metric, the set of spheres with common center x € R" is a 
totally geodesic subspace of E>i. The reason is that it is the fixed point set of a 
group of isometries acting on namely the group of rotations of R n around x. 
(We also have to assume uniqueness of solutions to the geodesic equation.) For the 
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G A metric and plane curves the set of concentric spheres has been studied in [16] 
and for Sobolev type metrics they have been studied in [5J |5] . Some work for the 
G°-metric has also been done by [15] . 

Theorem. Within a set of concentric spheres, any sphere is uniquely described by 
its radius r. Thus the geodesic equation within a set of concentric spheres reduces 
to an ordinary differential equation for the radius. It is given by: 



The space of concentric spheres is geodesically complete with respect to a G* metric 
iff 



and 



71— L I I I III * I , \ / \ 7 

r 2 Y $ ( r(i + g) '~^~ 1 ^ r J 7-1 



'" 1 ' ; l r(i + ny -( n - 1 )/ r ) dr = °° r 0>0- 



lr y lJ - V-"- ' 2 

For the metrics studied in this work, this yields: 

k to 

$ = Vol fc = — —r-r fe (™- 1 ) : incomplete 

r(l+f) fc 



Tt7T 2 i 



$ = 1 + A Tr(L) 2k = 1 + A {n jr : complete iff k > 

i±n .TrfL) 2 C(n) 
$ = Vol 1 -" +A ' = —~ : complete. 
Vol r Tl+1 



incomplete 

n + 1 



Proof. The differential equation for the radius can be read off the geodesic equation 
in section |6.2[ when it is taken into account that all functions are constant on each 
sphere, and that 

— n— 1 

Vol= r(l + f) ' i = -? Id ™' Tr(L fc ) = (-1)^. 

To determine whether the space of concentric spheres is complete, we calculate the 
length of a path / connecting a sphere with radius ro to a sphere with radius r\ : 



LfAf)= I ^/Gf{ffjf)dt 



$(Vol,Tr(L))r 2 vol(g)<it 

M 



/ kt W* ^ r ,-(n-l)/r)^ v dt 

J 1 "V vr(i + f)' v " /r(i + f) 
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r(i + 



— (n — l)/r jdr 



a 



10. Special cases of almost local metrics 
10.1. The G°-metric. The G°-metric is the special case of a G*-metric with $ 



1. Thus its geodesic equation can be read off from section 5.1 It reads as 



f t = a. V + Tf.fJ 

fu = ~\ (||/ t || 2 Tr(L> + T/.grad 9 (||/ t || 2 )) + (Tr(L).a - div" (hfj)) .f t . 



We have three conserved quantities, namely: 



5 (/ t T )voi( ff )er(T*M$ 


5 M vol(M)) 


reparametrization momentum 




/ /* vol(<7) 


linear momentum 


/ (/A/ t )vol(.g)eA 2 D 

J M 


r = so(n)* 


angular momentum 



The geodesic equation on Bi(M,R n ) is well studied. We can read it off from 
section [6j 



Tr(L).a 2 
1 ' 



ft = a.v, at = 
Sectional curvature is given by 

Ra(ai ) a 2 ,<i2,ai) = - / ||aida 2 - a 2 dai\\ 2 i vol(g) > 0. 
2 Jm 



This formula is in accordance with |15l section 4.5] since we have codimension one 
and a flat ambient space, so that only term(6) remains, and for the case of plain 
curves, it is in accordance with |17l section 3.5]. 

The G°-metric induces vanishing geodesic distance, see section [8j 



10.2. The G^-metric. For a constant A > 0, the G^-metric is defined as 

Gf(h,k)= f (l + ATr(Lf)g(h,k)vo\(g). 

This metric has been introduced by [THl HU HZ] • It corresponds to an almost local 
metric G* with §>(x,y,z) = (1 + Ay 2 ), thus its geodesic equation on Imm(M,R n ) 



ALMOST LOCAL METRICS ON SHAPE SPACE OF HYPERSURFACES IN n-SPACE 45 

is given by (see section 5.1 1: 



f t = a.u + Tf.fJ 1 



ftt — 



1 



A((2ATr(L))\\f t 
1 + ATt(L) 2 



2ATi(L 2 ) 



Tf. 



(2ATr(L))||/ t || 2 grad s (Tr(L))~grad 9 ((l + ATr(L) 2 )||/ t || 2 ) 



■ (^Tr(L)) 
.1 + ATr(L) 

+ div 9 (/ t T )-Tr(L). a 



2(l + ATr(L) 2 ) 
2 ; Aa + aTr(L 2 ) + dTr{L){fJ)) 

"ft- 



The conserved quantities have the form 



(1 + ylTr(L) 2 ) g{fj) vol(ej) € T(T*M ® M vol(M)) reparam. momentum 

/ (1 + ^4Tr(L) 2 )/ t vol(g) linear momentum 
Jm 

(1 + ATr(L) 2 )(f A ft) vol(.g) e /\ 2r " - «o(n)* angular momentum 



The horizontal geodesic equation for the G A metric reduces to 



ft = a-v 

a t = ia 2 Tr(L) 



2 AA(Tr(L)) +4Aag- 1 (dTr(L),da) 



(1 + ATt(L) 2 



2ATr(L)\\da\\ z g - 1 — ATv(L)Tv(L 2 )a 2 
(l + ATr(i) 2 ) 



For the case of curves immersed in R 2 , this formula specializes to the formula given 
in [ini section 4.2]. (When verifying this, remember that A = —D 2 in the notation 
of [E].) 

The curvature tensor Ro(ai, a2, oti, (12) is the sum of: 



P\P Q3Q3 negative terms, 



PqP Q5Q5 positive terms, and 
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Q1Q5 Q3Q5 indefinite terms. 



Ro(a%, a 2 , a%, a 2 ) — / A(ai Aa 2 - a 2 Aai) 2 vol(g) 
Jm 

+ / 2ATr{L)g2{{aida 2 — a 2 da\) ® {a\da 2 ~ a 2 da\), s) Yo\{g) 
Jm 



M 



1 

1 + ATi(L) 2 



— AA 2 g 1 (c?Tr(L) , a\da 2 — a 2 detij 



+ ATr(L) 2 ) 2 + 2A 2 Tr(L)A(Tr(£)) + 2A 2 Tr(L 2 ) Tr(L)'' 



\a\da 2 — a 2 dai\\ ! + (2A 2 Tr(X) 2 ) \\da± A da 2 || „2 



+ (8A 2 Tr(L))g^(dTr(L) ® (aida 2 - a 2 dai),dai A da 2 ) 



vol(g) 



We want to express the curvature in terms of the basic skew symmetric forms. 
Therefore, mimicking the notation of [Ttfl fT7] we define 

W2 = a\da 2 — a 2 da\, W 22 = a\Aa 2 — 0,2 Aai, W12 = dai A da 2 . 

Then the above equation reads as: 



Ro(ai, a 2l ai,a 2 ) = 


AW 2 2 2 vol(#)+ f 2ATr(L)^(W 2 ® W 2 , 


s) vol(s) 


JM 


Jm 


+ / M l + ATr(L)2 


-AA 2 g- x (dTi{L),W 2 ) 2 




-(l(l + ATr(L) 2 ) 


2 + 2A 2 Tr(L)A(Tr(L)) + 2A 2 Tr(L 2 ) Tr(L) 2 ) ||W 2 ||^_ 1 


+ (2A 2 Tr(L) 2 ) ||Wi 2 ||2g + (8A 2 Tr(L)) 5 °(dTr(L) <g> W 2 , H^ia) 


vo%) 



For the case of plain curves, this formula specializes to the formula given in |17l 
section 3.6]. 



The G^-metric satisfies condition (JlJ from section 8.6 
vanishing geodesic distance. 



thus it induces non- 



10.3. Conformal metrics. The conformal metrics correspond to almost local met- 
rics G* with $ = $(Vol). For the case of planar curves these metrics have been 
treated in [231 [Ml US Qj|] . [TH] provides very interesting estimates on geodesic dis- 
tance induced by metrics with $(Vol) = Vol and e Vo1 . The geodesic equation on 
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\mm{M, R") is given by: 



f t = h = a.v + Tf.h x , 



fc = - 



||/i|r vol(sf) Tr(L).i/ 



+ ||/ l || 2 Tr(L).^ + T/.grad 9 (||/ l || 2 ) 
$ v./.w 



Tr(i).avol(g)J + Tr(L).a - div 9 (/i T ) 



The conserved quantities are given by 



$(Vol)(;(/ t T ) vol(g) € T(T* M ®m vol(-M)) reparam. momentum 



$(Vol) / / t vol(<?) 



linear momentum 



3>(Vol) / (/A / t ) vol(g) € /\ 2 IR™ = so(n)* angular momentum 
7m 



The horizontal part of the geodesic equation is given by 
at = g(^H{a.u 1 a.v) — K(a.i/, a.v) 1 i^j 

= ~ (J « 2 vol( 5 )) Tr(L) + \a 2 Tr(L) + ^ (J a. Tr(L) vol(g) } a 

To simplify this equation let b{t) = <£>(Vol).a(t). We get 
b t = $'.(£>(/,„.„) Vol).a + $.a t 

= -$'.a. / Tr(L).a.vol( ff ) + $^a 2 .Tr(i) 

{ f a 2 vol( 5 ) ) . Tr(L) + $'.a. / Tr(L).a vol( 5 ) 
2 \JM J Jm 

= --$' / a 2 vol(.g).Tr(i) + -$a 2 .Tr(L). 
2 Jm 2 

Thus the geodesic equation of the conformal metric G* on B>i is 




For the case of curves immersed in R 2 , this formula specializes to the formula given 
in [13 section 3.7]. 

Assuming that $' and $" are non-negative, the curvature tensor consists of the 
following summands. 



Q2Q2 are the positive summands. 



P\P Q4Q4 Q1Q2 are the negative summands. 
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Q2Q4 is indefinite, but assuming that is a non-negative constant, it is negative. 

Solving the ODE ^ = C > leeds to $(Vol) = e C Vo1 . In the case of curves, 
conformal metrics of this type have been studied by [TJ] and |19j . 

is indefinite. 



Since the formula for sectional curvature with general $ = <£>(Vol) is still too 
long, we will only print the formula for <£>(Vol) = Vol. To shorten notation we will 
write a for the integral over a g C°°(M), i.e. 

a= / avol(g). 

J M 

Then the sectional curvature reads as: 



R (ai,a 2 ,ai,a 2 ) = -- Vol / \\aida 2 - a 2 da 1 \\ 2 , vo\(g) 
1 Jm 



4 Vol 
1 



a\.a 2 



a\ .TT(L) 2 al - 2aI^2~.Tr(L) 2 ai.a 2 + a 2 .Tr(L) 2 a 2 



4 Vol 
1 



(a 2 .Tr(L)a 2 2 



2ai.a 2 .Tr(L)ai.Tr(L)a 2 + a 2 .Tr(L)ai 



a 2 .Tr 9 ((da 2 ) 2 ) - 2aT^ .Tr 9 (d ai .da 2 ) + a|Tr 9 ((dai) 2 
- (af.a 2 .Tr(L 2 ) - 2.aYM^.ai.a 2 . Tr(L 2 ) + a|.a 2 .Tr(L 2 ) 



For the case of curves immersed in M 2 , this formula is in accordance with the 
formula given in [171 section 3.7]. 



From Condition ^ in section 8.6 we read off that the conformal metrics induce 
non- vanishing geodesic distance if $(Vol) > C. Vol for some constant C > 0. 



10.4. A scale invariant metric. For a constant A > we define the metric 

G s /(h,k) = J (Vol^ +A T ^f-)g(h,k)vo\(g). 

Scale invariance means that this metric does not change when /, h, k are replaced 
by A/, Xh, Xk for A > 0. To see that G is scale invariant, we calculate as in 
|17j how the scaling factor A changes the metric, volume form, volume and mean 
curvature. We fix an oriented chart (it , . . . , ii"^ 1 ) on M, Then 

(\frg(d i ,9 j )=9(T(^f)-9i,T(Xf).d j ) = X 2 .f*g(d l , dj) 

vo\((X.f)*g) = ^det(X 2 (f*g)\u) du 1 A ... A d^ 1 = A"" 1 vol( fg) 

Tr(L((A/)*s))) = ((A/)*5) ii 5(^^,^- / ) 

= 4(/*^(|^,^) = yTr(L(/)). 
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The scale invariance of the metric G SI follows. Thus along geodesies we have an 



additional conserved quantity (see section 5.2), namely: 



A I 



,i±n Tr(L) 2 



Vol 1 -" +A y - )g(f, ft) v ol(<?) scaling momentum 



From [6] we can read off the geodesic equation for G SI on Bi 
ft = a.v, 



a t = \a 2 Tr(L) 

Vol +A-^ 



2 v ; MSxi ^)' 



Vol 



l-n ^r„i2 



Tra) 2 \ f / •. „ N ,Tr(L 2 )Tr(L) 



For the case of curves immersed in K 2 , this formula specializes to the formula given 
in jTTJ, section 3.8]. (When verifying this, remember that A = — D 2 S in the notation 

of nao 

The metric G SI induces non-vanishing geodesic distance. This follows from the 
fact that log(Vol) is Lipschitz, see [111 section 3.8]. 

11. Numerical results 

11.1. Discretizing the horizontal path energy. We want to solve the boundary 
value problem for geodesies in shape space of surfaces in M 3 with respect to several 
almost local metrics, more specifically with respect to G*-metrics with 

$ = Vol fe , $ = e Vo1 , $ = 1 + A Tr(L) 2k 



1+3 . Tr(L) 2 



and the scale-invariant metric 

$ = VoIt^s +A 

Vol 

In order to solve this infinite-dimensional problem numerically, we will reduce it to 
a finite-dimensional problem by approximating an immersed surface by a triangular 
mesh. 

One approach to solve the boundary value problem is by the method of geodesic 
shooting. This method is based on iteratively solving the initial value problem for 
geodesies while suitably adapting the initial conditions. 

Another approach, and the approach we will follow, is to minimize horizontal 
path energy 

E hm \f)= [ [ $(Vol,Tr(i))g(/ t ,z,) 2 vol( ff ) 
Jo Jm 

over the set of paths / of immersions with fixed endpoints. Note that by definition, 
the horizontal path energy does not depend on reparametrizations of the surface. 
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Nevertheless we want the triangular mesh to stay regular. This can be achieved by 
adding a penalty functional to the horizontal path energy. 



11.2. Discrete path energy. To discretize the horizontal path energy 

B hor (/)= / f $(Vol,Tr(L)).g(/ t ,i,) 2 vol( 5 ), 
Jo Jm 

one has to find a discrete version of all the involved terms, notably the mean 
curvature. We will follow [20] to do this. Let V, E, F denote the vertices, edges and 
faces of the triangular mesh, and let star(p) be the set of faces surrounding vertex 
p. Then the discrete mean curvature at vertex p can be defined as 

1 1 vector mean curvature 1 1 | |V P (surface area) || 
' ||vector area|j || V p (enclosed volume)|| 

Here V p stands for a discrete gradient, and 



(vector mean curvature)^ = V p (surface area) = (cot ai + cot jii){p — pi) 

(p,Pi)eE 

is the vector mean curvature defined by the cotangent formula. In this formula, aj 
and fa are the angles opposite the edge (p,Pi) in the two adjacent triangles. For the 
numerical simulation it is advantageous to express this formula in terms of scalar 
and cross products instead of the cotangents. Furthermore, 

(vector area) p = V p (enclosed volume) = "(/)• (surface area of /) 

/eStar(p) 

is the vector area at vertex p. 
We discretize the time by 

= h < . . . < t N+1 = 1. 

Then the (N — l)(#V) free variables representing the path of immersions / are 

f(ti,p), with 2 < i < N, p G V. 

f(0,p) and f(l,p) are not free variables, since they define the fixed boundary 
shapes, ft can be approximated by either forward increments 



ffW(. s _ f(tj+l,P) ~ f(ti,p) 



ti+1 ~ H 
or backward increments 



tj t 



i-1 



We use a combination of both to make path energy symmetric. (Instead of this we 
could have also used the central difference quotient. However minimizing an energy 
functional depending on central differences favors oscillations, since they are not 
felt by the central differences.) Using the discrete definitions of normal vector and 
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increments we can calculate ff- at every vertex p and are now able to write down 
the discrete horizontal path energy: 



N t i 

H+l — U 



G*(^^) = ££$(vol,Tr(L)(p) 

p£V F3p 

.g(h(p),v(F)).g(k(p),v{F)) 
N t 

B te (/) = E ! 



area(stax^(p)) 
3 



i=l 



.(G* ti) (/r(tOj/ w (iO)+G* (ti+l) (/^fe +1 ),/ t ta (t 4+1 ))). 



This is not the only way to discretize the energy functional. There are several 
ways to distribute the discrete energy on faces, vertices and edges. Depending 
on how this was done, the minimizer converged faster, slower or even not at all. 
However if the minimizer converged to a smooth solution, the results were qualita- 
tively the same. This increased our belief in the discretization. However we do not 
guarantee the accuracy of the simulations in this section. 

This energy functional does not depend on the parametrization of the surface 
at each instant of time. So we are free to choose a suitable parametrization. Wc 
do this by adding to the energy functional a term penalizing irregular meshes. So 
instead of minimizing horizontal path energy, we minimize the sum of horizontal 
path energy and a penalty term. The penalty term measures the deviation of angles 
from the "perfect angle" 2n divided by the number of surrounding triangles, i.e. 



N 

EE E 

4=2 p£V (p,q,r)eA 



<(pq,pr) — (perfect angle) 



k e N. 



11.3. Numerical implementation. Discrete path energy depends on a very high 
number of real variables, namely three times the number of vertices times one 
less than the number of time steps. In the numerical experiments that we have 
done, this were between 5.000 and 50.000 variables. To solve this problem we 
used the nonlinear solver IPOPT (Interior Point OPTimizer [12]). IPOPT uses a 
filter based line search method to compute the minimum. In this process it needs 
the gradient and the Hessian of the energy. IPOPT was invoked by AMPL (A 
Modeling Language for Mathematical Programming [6] ). The advantage of using 
AMPL is that it is able to automatically calculate the gradient and Hessian. The 
user only has to write a model and data file for AMPL in a quite readable notation. 
The data file containing the definition of the combinatorics of the triangle mesh 
was automatically generated by the computer algebra system Mathematica. As an 
example, some discretizations of the sphere that we used can be seen in figure [I] 

11.4. Scaling a sphere. In section [9] we studied the set of concentric spheres in n 
dimensions. In dimension three the geodesic equation for the radius simplifies to 

r« - -r 2 t \ f 1 * + d^4r 2 7T + i(9 2 $) + \(d 3 <t>)] . 
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Figure 1. Triangulations of a sphere with 320, 500 and 720 tri- 
angles, respectively. 



This equation is in accordance with the numerical results obtained by minimizing 
the discrete path energy defined in section As will be seen, the numerics 

show that the shortest path connecting two concentric spheres in fact consists of 
spheres with the same center, and that the above differential equation is (at least 
qualitatively) satisfied. Furthermore, in our experiments the optimal paths obtained 
were independent of the initial path used as a starting value for the optimization. 

In all numerical experiments of this section we used 50 timesteps and a triangu- 
lation with 320 triangles. 

For conformal metrics of the type <I> = Vol fe and $ = e Vo1 , the differential 
equation for the radius is: 

$ = Vol fc 
$ = e Vol 

Note that the equation for $ = Vol -1 is r tt = 0. These equations have explicit 
analytic solutions given by 

r = C 1 ((k + 2)t-Ca) 1 ^' 

r=^y / log {C x t + C 2 ). 

A comparison of the numerical results with the exact analytic solutions can be 
seen in figure [2] and [3J The solid lines are the exact solutions. For the numerical 
solutions, 50 time steps and a triangulation with 320 triangles (see figure [I]) were 
used. Note that for big radii as in figure [5J the solution for 4> = e Vo1 has a very 
steep ascent, is more curved and lies above the solutions for $ = Vol, Vol 2 , Vol 3 . 
For small radii, it lies below these solutions, as can be seen in figure [3] Note also 
that when the ascent gets too steep, the discrete solution is somewhat inexact as 
in figure [2] 

For mean curvature weighted metrics, the differential equation for the radius is: 



rtt 
rtt 



,k + l 



\r 



Ave 



$ = Vol fc : 
$ = e Vo1 : 




Figure 3. Geodesies between concentric spheres of radius 0.1 to 
0.2 for several conformal metrics. Solid lines are the exact solu- 
tions. 
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Figure 5. Geodesies between concentric spheres for $ = 1 + 
ATt(L) 2 and varying A. Solid lines are the exact solutions. 
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r 

0.2 

0.18 



0.16 



0.14 



0.12 



0.2 



0.4 



0.6 



0.8 



$ = VoP 2 + A Tt(L) 2 Vol" 1 



Figure 6. Geodesies between concentric spheres for the scale- 
invariant metric. 



The numerics for these metrics are shown in figure [4] and figure [5j Note that we 
got convergence to a path consisting of concentric spheres even for the G°-metric 
(A = 0), even though we know from the theory that this is not the shortest path. 
In fact, there are no shortest paths for the G° metric since it has vanishing geodesic 
distance [T5] , 

For the scale-invariant metric, the differential equation is given by: 



* = Von 2 + ^ 
Vol 



Ttt 



<J_ 

r 



This equation has an explicit analytical solution: 

Tr(i) 2 _ 



$ = Vol" 



-A- 



Vol 



= de C2t . 



Note that this equation and therefore its solution is independent of A. Again, this 
is confirmed by the numerics, see figure [6] 



11.5. Translation of a sphere. In this section we will study geodesies between a 
sphere and a translated sphere for various almost local metrics of the type $ = Vol* , 
$ = e Vo1 and $ = 1 + ATr(L) 2k . 

Depending on the distance (relative to the radius) of the two translated spheres, 
different behaviors can be observed. 



High distance: 
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Shrink and grow: For some metrics it is possible to shrink a sphere in 
finite time to zero. For these metrics long translation goes via a shrinking 
and growing part. Metrics with this behavior are: $ = Vo\ k , $ = e Vo1 
and <f> = 1 + ATr(L) 2 . This phenomenon is studied in more detail in 
section [i 1.6 see also figure [8] 

Moving an optimal middle shape: For some of the metrics translation of 
a sphere with a certain optimal radius is a geodesic. For these metrics 
geodesies for long translations scale the sphere to the optimal radius and 
translate the sphere with the optimal radius. Metrics with this be havio r 
are $ = 1 + ATr(L) 2k for k > 1. This behavior is studied in section 11.7 



Low distance: 



Geodesies of pure translation. (<fr = 1 + A Tr(L) for k > 1, c.f. figure 11 ) 



Geodesies that pass through an ellipsoid, where the longer principal axis is 
in the direction of the translation (Conformal metrics, c.f. figure [7]). 
Geodesies that pass through an ellipsoid, where the principal axis in the 
direction of the translation is shorter ($ = 1 + j4Tr(L) 2fc for k > 1, c.f. 



figure 1 1 ) 



Geodesies that pass through a cigar shaped figure ($ — 1 + j4Tr(L) , c.f. 



figure 10 



Figure 7. Geodesic between two unit spheres translated by dis- 
tance 1.5 for $ = Vol. 20 timesteps and a triangulation with 500 
triangles were used. Time progresses from left to right. Boundary 
shapes t = and t = 1 arc not included. 



11.6. Shrink and grow. In section [9] we showed that it is possible to shrink a 
sphere to zero in finite time for some of the metrics, namely conformal metrics 
with $ = Vol fc or $ = e Vo1 and for the G A metric. For these metrics geodesies 
of long translation will go via a shrinking and growing part, and almost all of the 
translation will be done with the shrunken version of the shape. An example of 
such a geodesic can be seen in figure [8j 

We could not determine numerically whether a collapse of the sphere to a point 
occurs or not. But the more time steps were used, the smaller the ellipsoid in the 
middle turned out. Also, the energy of the geodesic path comes very close to the 
energy needed to shrink the sphere to a point and blow it up again. It is remarkable 
that almost all of the translation is concentrated at a single time step, independently 
of the number of timesteps that were used. The reason for this behavior is that 
high volumes are penalized so much: In the case of figure [sj e Vo1 is more than 1000 
times smaller in the middle than at the boundary shapes. 
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Figure 8. Geodesic between two unit spheres translated by dis- 
tance 2 for <£> = e Vo1 . 20 timesteps and a triangulation with 500 
triangles were used. Time progresses from left to right. Boundary 
shapes t = and t — 1 arc not included. 



We now want find out under what conditions on the distance and radius of the 
boundary spheres of the geodesic this behavior can occur. To do this, we compare 
the energy needed for a pure translation with the energy needed to first shrink the 
sphere to almost zero, then move it, and then blow it up again. 

The energy needed for a pure translation of a sphere with radius r by distance 
I in the direction of a unit vector e\ is given by 

E= [ [ <f>(Vol,Tr(L))g(e.ei,v) 2 vol(g)dt 
Jo Js 2 



2 /* 7T r2TT I I COS tp Sill \j i i 

— $(4r 2 7r, — ) J J g \i.ei, I sintpsinfl I I r 2 sin 9 dtpdO 




2 



/0 Jo \ \ COS0 

= $(4r 2 7r, — )/ / (cos </> sin 60 V sin 6^(2(9 = $(4r 2 7r, — ).— ^ 2 .r 2 
r Jo Jo ' r 3 

Any other unit vector can be chosen instead of ei, yielding the same result. 




■ $ = Vol $ = Vol 2 $ = Vol 3 ■ $ = e Vo1 ■ G 1 



Figure 9. Left: Shrinking a sphere to zero along a geodesic path 
and blowing it up again. Right: Pairs of £ and r such that trans- 
lating a sphere of radius r by distance t needs as much energy as 
shrinking it to zero and blowing it up again. G 1 stands for the G A 
metric with A = 1. 



We will now calculate the energy needed for shrinking the sphere, moving it, and 
blowing it up again. The energy needed for translating a sphere of radius almost 
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zero can be neglected. Shrinking and blowing up is done using the solutions to the 
geodesic equation for the radius from the last section, where one has to adapt the 
constants to the boundary conditions. For the shrinking part, we have r(0) = r and 
r(g) =0, and for the growing part we have r(|) = 0, r(l) = r, see figure M (left) . 

The energy of the path is 

r 1 r Ak+2 k+i 
$ = Vol fc : E = / Vol fc / r 2 vol(g)dt = — ^r 2fe+4 



/o Js 2 (k 

$ = e Vo1 : E= C e Vo1 / r\ vol(g)dt = -(e 2 ^ 2 - l) 2 . 

Jo Js 2 7r 

The energy of the two different paths are the same when 

* = Vol*: i^-fl- 

k + 2 

$ = e v ol; i= V3(l-e~^)_ 

2rir 

These curves are shown in figure [9] (right). We did not derive an analytic solution 
for the G A metric, but for A = 1 one can see the solution curves in figure [9] 

11.7. Moving an optimal shape. In the following we want to determine whether 
pure translation of a sphere is a geodesic. Therefore let ft — fo + b(t) ■ ei, where 
/o is a sphere of radius r and where b(t) is constant on M. Plugging this into the 



geodesic equation from section 5.1 yields an ODE for b(t) and a part which has to 



vanish identically. The latter is given by: 

(1) (di$)-4rV + (9 2 $)^+$- =0 

For conformal metrics this equation is only satisfied if $ = Vol - . Since this metric 
induces vanishing geodesic distance (see section [8]) we are not interested in this 
case. For curvature weighted metrics the above equation reads as: 

$ = l + ATr(L) 2fc : ~ = 1 

Solutions to these equations are given by: 

^ = 1 + ATr(L) 2k : r = 2 2 {/A(k-l), k > 1. 

For the most prominent example the G A metric this yields r = and therefore 
translation can never be a geodesic for this type of metrics. The numerics have 
shown that the G A metrics yields geodesies that resemble the geodesies of the G A 
metric for planar curves from |16l section 5.2]. Namely, when the two spheres 
are sufficiently far apart, the geodesic passes through a cigar-like middle shape, 



see figure 10 As predicted by the theory (see section 8.8) geodesies for very high 



distances tend to have a similar behavior as Vol fe metrics, i.e. the geodesic first 



shrinks the sphere, then moves it, and then blows it up again (cf. section 11.6 1 



For metrics weighted by higher factors of mean curvature the above equation 
for the radius has a positive solution. For these metrics geodesies for translations 
tend to scale the sphere until it has reached the optimal radius and then translate 
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Figure 10. Middle figure of a geodesic between two unit spheres 
translated by distance 3 for $ = 1 + ATr(L) 2 . From left to right: 
A = 0.2, A = 0.4, A = 0.6, A = 0.8. In each of the simulations 20 
timesteps and a triangulation with 720 triangles were used. 

it. If the radius is already optimal the resulting geodesic is a pure translation (see 
figure 11). 
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Figure 11. Geodesic between two unit spheres translated by dis- 
tance 3 for $ = 1 + Jg Tr(L) 4 (first row) and $ = l+Tr(L) 6 (second 
row) . In each of the experiments 20 timesteps and a triangulation 
with 720 triangles were used. Time progresses from left to right. 
Boundary shapes t — and t — 1 are not included. 



If the distance is not high enough there still occurs a scaling towards the optimal 
size, but the middle figure is not a perfect sphere anymore. Instead it is an ellipsoid 
as in figure |11| 



Figure 12. Geodesic between a sphere and a sphere with a small 
bump for $ = Vol. 20 timesteps and a triangulation with 500 
triangles were used. Time progresses from left to right. 



(ill 



MARTIN BAUER, PHILIPP HARMS, PETER W. MICHOR 



11.8. Deformation of a shape. We will calculate numerically the geodesic be- 
tween a shape and a deformation of the shape for various almost local metrics. 
Small deformations are handled well by all metrics, and they all yield similar re- 
sults. An example of a geodesic resulting in a small deformation can be seen in 



figure 12, where a small bump is grown out of a sphere. The energy needed for this 
deformation is reasonable compared to the energy needed for a pure translation. 
Taking the metric with $ = Vol as an example, growing a bump of size 0.4 as in 
figure fl2] costs about a third of a translation of the sphere by 0.4. 

Bigger deformations work well with Vol fe -metrics and curvature weighted metrics, 
but not with the e Vol -metric, which tends to shrink the object and to concentrate 
almost all of the deformation at a single time step. In figure [13} a large deformation 
can be seen for the case of <& = Vol and $ = e Vo1 . Clearly one can see that the 
e Vol -metric concentrates almost all of the deformation in a single time step. We 
have met this misbehavior of the e Vol -metric already with translations. Again, the 
reason is that e Vo1 is so sensitive to changes in volume. 







^^^^^^^ 



C 



Figure 13. Large deformation of a shape for <f> = Vol and $ = 
e Vo1 . 20 timesteps and a triangulation with 500 triangles were used. 
Time progresses from left to right. 



In figure [14| one sees that higher curvature weights smoothen the geodesic. 
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FIGURE 14. Large deformation of a shape for $ = 1 + 0.1 Tr(L) 2 
(top), $ = l + 10Tr(i) 2 (middle) and $ = l + Tr(L) 6 (bottom). 20 
timesteps and a triangulation with 720 triangles were used. Time 
progresses from left to right. 
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The AMPL model file 



Listing 1. AMPL model file 

l param A default 1; 

param k default 1; 
3 param B default 1; 

param 1 default 1; 
5 param TimestepsN > 1 integer; 

param VerticesN integer; 
7 param Penalty Factor default 1; 

param PenaltyExponent default 2; 
9 set VerticesI := 1.. VerticesN; 

set VerticesOfEdgesI within {VerticesI , VerticesI }; 
ll set VerticesOfFacesI within { VerticesI , VerticesI , VerticesI } ; 

set FacesOfVerticesI {v in VerticesI} within VerticesOfFacesI; 
13 set LinkOfVerticesI {VerticesI} within {VerticesOfFacesI, VerticesOfEdgesI, { — 1,1}}; 

set AdjacentEdgesOfVerticesI {VerticesI} within {VerticesOfEdgesI, {1,-1}, VerticesOfEdgesI, {1,-1 
15 set EdgesOfFacesI {VerticesOfFacesI} within VerticesOfEdgesI; 

set EdgesOfVerticesI {v in VerticesI} := setof {(fl , f2 , f3 ,el,e2,o) in Link0fVerticesI[v]}(el,e2); 

17 

param Pi default 3.141592653589793; 
19 param PerfectAngle {v in VerticesI} default cos(2*Pi/card(FacesOfVerticesI[v])); 

param InitialVertices {VerticesI ,1..3}; 
21 param FinalVertices {VerticesI, 1.. 3}; 

23 var Middle Vertices {2..TimestepsN,VerticesI ,1..3}; 

25 var Vertices {t in 1.. TimestepsN+l,v in VerticesI, i in 1..3} = 

(if t=l then InitialVertices [v, i] 
27 else if t=TimestepsN+l then FinalVertices[v,i] 



else Middle Vertices [t ,v, i ]); 

29 

var VectorOfEdges {t in l..TimestepsN+l, (vl,v2) in VerticesOfEdgesI,i in 1..3} = 
31 Vertices [t,v2,i] — Vertices[t ,vl,i ]; 

33 var LengthOfEdges {t in l..TimestepsN+l, (vl,v2) in VerticesOfEdgesI} = 

sqrt(VectorO£Edges[t,vl,v2,l]"2+VectorOfEdges[t,vl,v2,2]"2+VectorOfEdges[t,vl,v2,3]"2); 

35 

var CrossOfFaces {t in l..TimestepsN+l,(vl,v2,v3) in VerticesOfFacesI, i in 1..3} = 
37 if i=l then (Vertices[t,v2,2]-Vertices[t,vl,2])*( Vertices [t,v3,3]-Vertices[t,vl, 3]) - 

(Vertices [t, v2,3]-Vertices[t ,vl,3])*( Vertices [t,v3,2]-Vertices[t,vl, 2]) 
39 else if i=2 then — (Vertices[t,v2,l]— Vertices[t,vl ,1])*( Vertices [t,v3, 3]— Vertices [t,vl, 3]) + 

(Vertices [t ,v2,3]— Vertices [t ,vl ,3])*( Vertices [t ,v3,l]— Vertices [t ,vl ,1]) 
41 else (Vertices [t, v2,l]-Vertices[t,vl,l])*( Vertices [t,v3,2]-Vertices[t,vl, 2]) - 
(Vertices [t, v2,2]-Vertices[t,vl,2])*( Vertices [t,v3,l]-Vertices[t,vl,l]) ; 

43 

var NormCrossOfFaces {t in l..TimestepsN+l,(vl,v2,v3) in VerticesOfFacesI} = 
45 sqrt(CrossOfFaces[t,vl,v2,v3,l]~2 + CrossOfFaces[t,vl,v2,v3,2]~2 + CrossOfFaces[t,vl,v2,v3,3] "2); 

47 var NuOfFaces {t in l..TimestepsN+l,(vl,v2,v3) in VerticesOfFacesI, i in 1..3} = 
CrossOfFaces[t,vl,v2,v3,i]/NormCrossOfFaces[t,vl,v2,v3]; 

49 

var AreaOfFaces {t in l..TimestepsN+l,(vl,v2,v3) in VerticesOfFacesI} = 
5 1 NormCrossOfFaces [t , v 1 , v2 , v3] /2 ; 

53 var AreaOfVertices {t in 1.. TimestepsN+1, v in VerticesI} = 

(sum {(fl,f2,f3) in FacesOfVerticesI[v]} AreaOfFaces [t,fl,f2 , f3 ])/3; 

55 

var VectorAreaOfVertices {t in l..TimestepsN+l, v in VerticesI, i in 1..3} = 
57 (sum {(vl,v2,v3) in FacesOfVerticesI[v]} CrossOfFaces[t,vl,v2,v3,i ])/6; 



59 var SquareOfNormOfVectorAreaOfVertices {t in l..TimestepsN+l, v in VerticesI} = 

VectorAreaOfVertices[t,v,l] ~2+VectorAreaOfVertices[t,v,2] ~2+VectorAreaOfVertices[t,v,3] "2; 

61 

var NormOfVectorAreaOfVertices {t in l..TimestepsN+l, v in VerticesI} = 
63 sqrt( SquareOfNormOfVectorAreaOfVertices[t,v]); 

65 var Volume {t in l..TimestepsN+l} = 

sum{(vl,v2,v3) in VerticesOfFacesI} AreaOfFaces[t,vl,v2,v3]; 

67 

var VectorMeanCurvatureOfVertices {t in l..TimestepsN+l, v in VerticesI, i in 1..3} = 
69 if i = l then 

sum {(fl,f2,f3,el,e2,o) in LinkOfVerticesI[v]} o* 
71 ( VectorOfEdges[t,el,e2,2]*NuOfFaces[t,fl,f2,f3,3] - 

VectorO£Edges[t,el,e2,3]*NuOfEaces[t,fl,f2,f3 ,2] ) 
73 else if i=2 then 

sum {(fl,f2,f3,el,e2,o) in LinkOfVerticcsI[v]} o* 
75 (-VectorOfEdges[t,el,e2,l]*NuO£Faces[t,fl,f2,f3,3] + 

VectorOfEdges[t,el,e2,3]*NuOfFaces[t,fl,f2,f3 ,1] ) 

77 else 

sum {(fl,f2,f3,el,e2,o) in LinkOfVerticesI[v]} o* 
79 ( VectorOfEdges[t,el,e2,l]*NuOfFaces[t,fl,f2,f3,2] - 

VectorOfEdges[t,el,e2,2]*NuOfFaces[t,fl,f2,f3,l] ) ; 

81 

var SquarcOfScalarMeanCurvatureOfVertices {t in l..TimestepsN+l, v in VerticesI} = 
83 (VectorMeanCurvatureOfVertices[t,v,l] "2+VectorMeanCurvatureOfVertices[t,v,2] "2 

+VectorMeanCurvatureOfVertices[t,v,3] "2)/SquareOfNormOfVectorAreaOfVertices[t,v] ; 

85 

var PhiOfVertices {t in 1.. TimestepsN+l,v in VerticesI} = 
87 1+ A*(SquareOfScalarMeanCurvatureOfVertices[t,v])~k +B*(Volume[t])"l; 



89 var IncrementsOfVertices {t in 1.. TimestepsN,v in VerticesI, i in 1..3} = 
TimestepsN*(Vertices[t+l,v,i] — Vertices[t , v, i ]); 

91 

var Energy = 1/ 12 / TimestepsN *( 
93 sum {t in l..TimestepsN,v in VerticesI} 

PhiOfVertices[t , v] * sum {(wl,w2,w3) in FacesOfVerticesI[v]} 5 
95 ( IncrementsOfVertices[t,v,l]*CrossOfFaces[t,wl,w2,w3,l] + ™ 

IncrementsOfVertices[t,v,2]*CrossOfFaces[t,wl,w2,w3,2] + 2 
97 IncrementsOfVertices[t,v,3]*CrossOfFaces[t,wl,w2,w3,3] )"2 / > 

NormCrossO£Faces[t,wl,w2,w3] + B 
99 sum {t in l..TimestepsN,v in VerticesI} 

PhiOfVertices[t+l,v] * sum {(wl,w2,w3) in FacesOfVerticesI[v]} S 
101 ( IncrementsOfVertices[t,v,l]*CrossOfFaces[t+l,wl,w2,w3,l] + 

IncrementsOfVertices[t,v,2]*CrossOfFaces[t+l,wl,w2,w3,2] + ^ 
103 IncrementsOfVertices[t,v,3]*CrossOfFaces[t+l,wl,w2,w3,3] )~2 / > 

NormCrossO£Faces[t+l,wl,w2,w3] ); | 
105 ^ 
var Penalty = H 
107 sum {t in l..TimestepsN+l, v in VerticesI, (vl,wl,ol,v2,w2,o2) in AdjacentEdgesOfVerticesI[v]} fa 
abs( * 
109 ( VectorOfEdges[t,vl,wl,l]*VectorO£Edges[t,v2,w2,l] + g 

VectorO£Edges[t,vl,wl,2]*VectorO£Edges[t,v2,w2,2] + § 
ill VectorOfEdges[t,vl,wl,3]*VectorOfEdges[t,v2,w2,3] ) * ol * o2 § 

/ LengthOfEdges[t,vl,wl] / LengthOfEdges[t,v2,w2] 
113 - Perfect Angle [v] 

) " Penalty Exponent ; 

115 

minimize f: 

117 Energy+Penalty*PenaltyFactor; 
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